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Summary  of  Technical  Effort 

This  report  summarizes  the  theoretical  development  of  a  nev 
approach  to  calculating  molecular  ground  and  excited  adiabatic  potential 
energy  surfaces  and  their  nonadiabatic  coupling  matrix  elements.  The 
fundamentals  of  the  method  will  be  briefly  described  herein;  the  details  can 
be  found  in  the  published  papers  resulting  from  our  investigations,  which  are 
included  as  Appendices  to  this  report. 

Methods  based  on  single-reference  many-body  perturbation  theory 
(MBPT)  have  been  shown  to  give  highly  accurate  results  for  a  large  number  of 
molecules.  These  methods  have  become  widely  used,  and  the  diagrammatic 
theory  has  provided  an  understanding  of  correlation  effects  that  has  led  to 
the  now  almost  universal  practice  of  adding  approximately  correcting  terms  to 
energies  calculated  by  configuration  interaction  methods.  However,  the 
single  model-function  methods  are  not  appropriate  for  many  systems  of 
chemical  interest  or  for  excited  states  in  general.  These  methods  cannot 
describe  the  breaking  of  multiple  bonds,  and  they  can  obtain  only  a  single 
state  of  a  given  symmetry.  Since  most  geometrical  arrangements  of  a 
polyatomic  system  do  not  have  any  symmetry,  this  means  that  no  excited  states 
may  be  calculated.  Furthermore,  the  proper  description  of  ground  states  of 
open  shell  species  often  requires  more  than  a  single  model  function.  A 
method  with  the  capability  of  treating  open  shell  systems  and  excited  states 
is  a  necessary  tool  for  the  theoretical  investigation  of  many  chemical 
phenomena.  In  turn,  the  basic  and  detailed  understanding  of  these  processes 
is  critical  to  many  areas  of  defense  technology. 

The  development  of  multidimensional  many-body  perturbation  theory 
retains  the  beneficial  characteristics  of  MBPT  and  extends  the  potential 
applicability  to  all  of  the  cases  described  above.  This  generalization 
ultimately  requires  a  prescription  for  the  energy  denominators  in  the 
diagrammatic  expansion.  In  place  of  the  single  zeroth-order  model-state 
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energy  required  in  MBPT,  a  set  of  differing  energies,  each  appropriate  for  an 
excited  state  of  interest,  must  be  incorporated.  There  are  several 
alternatives  for  this.  We  have  chosen  to  base  our  method  on  a  canonical  van 
Vleck  formalism.  This  leads  to  equal  treatment  of  the  states  of  interest, 
which  are  obtained  from  a  single  diagonal ization  of  an  inherently  symmetric 
effective  Hamiltonian  matrix.  This  formalism  is  presented  in  detail  and 
compared  with  other  approaches  in  Appendix  A.  Our  efforts  in  this  phase  of 
development  included  the  adoption  of  a  concise  notation  which  simplified  the 
description  of  our  formalism,  as  well  as  that  of  others.  This  was  important 
because  so  many  different  notations  and  diagrams  have  been  employed  by 
various  authors  that  the  similarities  and  relationships  among  their 
approaches  have  often  been  obscured. 

We  chose  to  implement  the  van  Vleck  formalism  analogously  to  the 
MBPT  method;  i.e.,  in  terms  of  diagrams  drawn  with  respect  to  an  SCF  (self- 
consistent  field)  determinant.  These  diagrams  provide  a  compact  visual 
representation  and  classif ication  of  the  terms  in  the  arithmetic  expressions 
for  each  order  of  the  effective  Hamiltonian  matrix  elements.  Diagrammatic 
rule  can  be  set  forth  which  make  them  entirely  sufficient  for  construction  of 
these  matrix  elements  to  any  order  for  an  arbitrary  set  of  model  space 
determinants.  Appendix  B  demonstrates  the  diagrammatic  implementation  of  the 
method  through  second  order  for  a  model  space  of  non-,  singly-,  and  doubly- 
excited  model  states. 

Traditionally,  quantum  chemists  and  chemical  physicists  could  be 
divided  into  two  mutually  exclusive  groups:  electronic  structure  specialists 
and  dynamicists.  So  far,  the  interaction  between  these  two  factions  has  been 
limited  to  the  former  providing  potential  energy  surfaces  (PES),  usually  for 
ground  states,  to  the  latter.  However,  many  chemical  phenomena  involve  more 
than  one  electronic  state.  Some  of  the  most  interesting  processes  involve 
the  avoided  crossing  of  states  of  the  same  symmetry.  The  accurate  simulation 


of  such  phenomena  requires  not  only  the  PES  but  also  the  nonadiabatic 
coupling  matrix  elements.  These  quantities  are  not  ordinarily  calculated; 
even  when  they  are,  a  problem  arises  from  their  coordinate  system  dependence 
because  the  structure  calculations  use  a  different  system  from  the  dynamical 
one.  The  last  phase  of  our  method  development  was  devoted  to  narrowing  this 
information  gap. 

Appendix  C  gives  the  details  of  a  novel  scheme:  the  diagrammatic 
perturbative  determination  of  nonadiabatic  coupling  matrix  elements.  This 
approach  prescribes  the  couplings  directly  in  terms  of  expressions  (or 
diagrams)  whose  components  consist  of  (1)  components  of  the  effective 
Hamiltonian  matrix  elements  (i.e.,  two-electron  integrals  and  orbital  energy 
differences)  and  (2)  a  set  of  one-electron-type  integrals.  With  this  method, 
the  wavefunction  is  never  explicitly  calculated  or  differentiated.  The 
perturbative  approach  prescribes  the  direct  construction  of  increasingly 
accurate  values  of  the  couplings.  Furthermore,  the  same  set  of  diagrams  can 
be  used  to  obtain  all  of  the  information  required  to  transform  between 
coordinate  systems  by  simply  substituting  different  one-electron  integrals. 
Thus,  the  theoretical  foundation  necessary  for  the  accurate  study  of 
dynamical  systems  using  ab  initio  electronic  states  and  couplings  has  been 
laid. 

In  addition  to  the  subjects  of  the  published  works  in  the 
Appendices,  considerable  effort  has  been  devoted  to  detailed  examination  of 
the  cancellation  of  unlinked  diagrams  and  the  form  of  the  renormalization 
terms  in  third  and  higher  orders.  We  have  also  begun  implementation  of  the 
method  into  a  computer  code.  Although  we  do  not  have  final  results  yet,  the 
results  of  similar  effective  Hamiltonian  methods  indicate  their  potential 
successful  application  to  systems  of  technological  interest.  To  our 
knowledge,  the  diagrammatic  perturbative  calculation  of  nonadiabatic 
couplings  is  unprecedented. 
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Three  forms  of  quasidegenerate  perturbation  theory  are  discussed  and  compared  in  terms  of  a  common 
general  formulation  based  on  a  similarity  transformation  which  decouples  the  model  space  and 
complementary  space  components  of  the  Hamiltonian.  The  discussion  is  limited  to  formal,  rather  than  many- 
body  (diagrammatic),  aspects.  Particular  attention  is  focused  on  a  "canonical"  form  of  van  Vleck 
perturbation  theory,  for  which  new  and  highly  compact  formulas  are  obtained.  Detailed  comparisons  are 
made  with  the  Kirtman-Certain-Hirschfelder  form  of  the  van  Vleck  approach  and  with  the  approach  based 
on  intermediate  normalization  which  has  been  used  as  the  basis  for  most  of  the  diagrammatic  formulations  of 
quasidegenerate  perturbation  theory. 
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canonical  van 


I.  INTRODUCTION 

Quasidegenerate  perturbation  theory  (QDPT)  has  been 
receiving  increasing  attention  in  recent  years  (for  sev¬ 
eral  excellent  discussions  and  reviews  see  Refs.  1-10). 
It  provides  the  perturbation  theory  analog  of  '‘multiref¬ 
erence”  configuration  interaction  (Cl)  techniques, 11,12 
which  have  proved  effective  in  the  treatment  of  potential 
surfaces  and  excited  states  of  molecular  systems. 
“One-dimensional"  perturbation  theory  (based  on  a  sin¬ 
gle-configuration  zero-order  function)  has  been  quite 
successful  in  applications  to  many  near-equilibrium 
ground  state  systems  and  in  some  other  cases  in  which 
a  single  configuration  provides  a  reasonably  adequate 
starting  point  (for  some  examples  see  Refs.  13-21). 
Practically  all  such  applications  have  used  the  many- 
body  form  of  Rayleigh-Schrodinger  perturbation  theory 
(RSPT)  (for  reviews  see,  for  example.  Refs.  22  and 
23).  This  approach  has  significant  advantages  over  Cl 
techniques  because  of  its  use  of  the  linked  cluster  ex¬ 
pansion24' 25  and  its  “extensivity”  property, 2,18,18,28,27 
i.  e. ,  its  correct  scaling  with  the  size  of  the  system. 
While  coupled  cluster  techniques2*  appear  capable  of 
extending  the  range  of  usefulness  of  single-configura¬ 
tion  based  treatments  considerably,21’2*  including  some 
nearly  degenerate  or  even  fully  degenerate  cases,30121 
it  still  appears  desirable  to  have  practical  computa¬ 
tional  techniques  based  on  a  multidimensional  “model 
space”  as  the  zero-order  approximation.  (An  alternative 
strategy  is  to  use  unrestricted  Hartree-Fock  zero-order 
functions, 1*’1*  but  this  appears  to  have  some  serious  dis¬ 
advantages.21)  In  fact,  QDPT  may  be  expected  to  pro¬ 
vide  faster  convergence  and  more  general  applicability 
than  one-dimensional  perturbation  expansions,  and 
should  become  an  increasingly  important  tool  for  the 
calculation  of  highly  correlated  electronic  wave  func¬ 
tions  and  energies  of  atoms  and  molecules.  (A  coupled 
cluster  analog  of  QDPT  has  also  been  formulated.1’32) 

In  the  present  contribution  we  are  concerned  primari¬ 
ly  with  a  particular  form  of  van  Vleck  perturbation  the- 
ory4'*’33-42  (VVPT)  which,  in  the  spirit  of  Klein’s  treat¬ 
ment*  (see  also  Jdrgensen42  and  Brandow4),  we  shall 


call  “canonical”  VVPT.  (Klein,  as  well  as  Primas, 38,38 
treated  the  case  of  exactly  degenerate,  rather  than 
quasidegenerate,  zero-order  subspaces. )  Combining 
some  elements  from  the  treatments  of  Primas38,39  and 
of  Jdrgensen, 40,41  we  shall  present  a  simple  derivation 
of  the  canonical  VVPT  formalism,  and  obtain  highly 
compact  expressions  for  the  decoupling  operator  and 
the  resulting  effective  Hamiltonian.  We  shall  compare 
this  formalism  with  another  version  of  VVPT  discussed 
by  Kirtman45,44  and  Certain  and  Hirschfelder,  45,48  and 
shall  present  simple  derivations  of  explicit  equations 
connecting  canonical  VVPT  with  the  more  common 
QDPT  formulation  based  on  intermediate  normaliza¬ 
tion  and  a  non-Hermitian  model  Hamiltonian. 2-10  The 
different  quasidegenerate  approaches  will  be  discussed 
in  terms  of  a  common  general  formulation  (see  also 
Klein5  and  Brandow4)  which  clearly  shows  their  relation¬ 
ships.  Only  the  formal  aspects  of  the  QDPT  expansions 
will  be  discussed  here.  We  expect  to  discuss  their 
many-body  realization  in  terms  of  diagrammatic  ex¬ 
pansions  in  future  contributions. 

The  notation  to  be  used  and  the  common  framework 
for  the  treatments  of  the  different  formalisms  are  pre¬ 
sented  in  Sec.  II.  Canonical  VVPT  is  derived  in  Sec. 

III.  The  Kirtman-Certain-Hirschfelder  (KCH)  form  of 
VVPT  and  the  intermediate  normalization  form  of  QDPT 
are  discussed  in  Secs.  IV  and  V,  respectively.  The  re¬ 
sults  are  discussed  in  Sec.  VI. 

II.  COMMON  FRAMEWORK  FOR  QUASIDEGENERATE 
PERTURBATION  FORMALISMS 

A.  Notation 

The  Hamiltonian  H  is  partitioned  into  a  zero-order 
part  and  a  perturbation 

H  =W0+  V  .  (1) 

The  eigenfunctions  of  f/0  will  be  written  in  the  form  If), 
with  eigenvalues  t,,  as 

#oif>  =*,!<>  .  (2) 
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and  the  set  of  these  eigensolutions  will  be  partitioned 
into  two  subsets 

{f }  =  l«,  *3, . . . }  U  {ij, . . . }  (3) 

defining  the  ‘‘model  space”  {a,  (3. ... }  and  its  orthogonal 
complement.  The  projection  operator  into  the  model 
space  is 

/>=£  |o)<«|  ,  (4) 

01 

and  its  orthogonal  complement  is 

y.i-p.£l0<f|  .  (5) 

i 

Any  operator  A  can  be  partitioned4’  into  a  block  diag¬ 
onal  part  Ad  and  a  block  off-diagonal  part  Ax 

A  -Ad  -  Ax  ,  (6) 

A0  =  PAP  +  QAQ  ,  Ax  =PAQ  +  QAP  .  (7) 

For  a  product  of  two  operators  we  have 

[AB)q  =Aq  B D  +  Ax  Bx  ,  (AB)X  ~Ap  Bx  +AX  Bq  .  (8) 


if„=X>l/3>Cs<,=U|a)  ,  (14) 

* 

where 

i<J>  =  £l/3  >C#a  (15) 

a 

are  the  “bonne  functions”  referred  to  by  Bloch48  and  by 
Jdrgensen  and  Pedersen. 40,41 

Rewriting  Eq.  (11)  in  the  form 

UK  =  HU  ,  (16) 

and  splitting  it  into  diagonal  and  off-diagonal  blocks  us¬ 
ing  Eqs.  (8)-(10),  we  find  that  the  condition  (12)  leads 
to  the  following  implicit  equations: 

HaUx=-VxUD+UxX  ,  (17) 

or,  using  the  definition  of  IF  in  Eq.  (12), 

[Ha,  Ux)=-VxUd-Vd  Ux  -  Vx  IF  ,  (18) 

and 

W=  [«„,  UD\  +  Vx  Ux  *  VD  VD  -  (C0  -  DIF  .  (19) 


For  the  Hamiltonian  we  note  that,  since  H0  is  diagonal, 
HD=H0+VD  ,  (9) 

HX  =  VX  .  (10) 

B.  Quasidegenerate  perturbation  theory 

The  essential  feature  of  the  various  QDPT  formalisms 
is  a  similarity  transformation  which  block  diagonalizes 
the  Hamiltonian 

3C  =  U-lHU  ,  (11) 

with 

K=X0=H0+W  ,  Kx=0.  (12) 

This  is  not  always  the  form  in  which  the  formalisms  are 
presented  (particularly  for  the  intermediately  normal¬ 
ized  form1-10),  but  it  provides  a  common  and  simple 
starting  point  for  straightforward  derivation  of  the  equa¬ 
tions  for  all  of  them.  The  decoupling  operator  U  is 
unitary  in  the  van  Vleck  formalisms  and  produces  a 
Hermitian  effective  (or  “model")  Hamiltonian  PKP.  It 
is  nonunitary  in  the  intermediate  normalization  approach 
(where  it  is  referred  to  as  the  "wave  operator")  and 
leads  to  a  non-Hermitian  jC  (which  can  easily  be  trans¬ 
formed  to  a  Hermitian  form,  if  desired4).  The  operator 
IF  defined  in  Eq.  (12)  is  referred  to  as  the  "level  shift” 
operator,  particularly  in  the  exactly  degenerate  case. 

Obviously,  X  has  the  same  eigenvalues  as  H,  so  that 
diagonalization  of  the  model  Hamiltonian  PKP  provides 
a  subset  of  the  eigenvalues  of  H.  The  perturbed  model 
functions  are 

=  ,  (13) 

( 

and  the  overlap  matrix  between  them  is  Pt/’t/P  (which  is 
a  unit  matrix  in  the  van  Vleck  case).  The  correspond¬ 
ing  eigenfunctions  t’a  of  H  are  obtained  by  transforming 
the  perturbed  model  functions  (13)  with  the  matrix  C  of 
eigenvectors  (right  eigenvectors  in  the  non-Hermitian 
case)  of  PKP: 


Expanding  the  operators  in  orders  of  the  perturbation 


O-Z 

tf<n\  x*£; 

(20) 

n»0 

f!»0 

n*  l 

with 

U(0>  =  1 

1p(0)  _  rr 

,  jv  -n0  , 

•JCfa>  =  IF*’”  (n>0)  , 

(21) 

explicit  recursive  equations  are  obtained  for  >i>  0: 

PI-1 

[»o,  tfr”)  =  -  Vx  U'S’"  -  VD  U?-v  +  £  U?"  IF0’-"”  , 


IF0”  =  [H„,  Ulo')*Vx  U?'l)  -  VD  UtT"  -  Y2  U(gm>  IV'4’1*”' 4  . 

1 

(23) 

However,  the  decoupling  operator  U  is  not  fully  deter¬ 
mined  by  the  condition  (12)  or,  equivalently,  by  Eqs. 

(22)  and  (23)  (note  that  no  equation  for  I'i”’  has  been  ob¬ 
tained).  In  fact,  multiplication  of  tT  by  any  block  diag¬ 
onal  operator  does  not  destroy  the  decoupling  of  X. 
Different  supplementary  conditions  on  V  (i.e. ,  specifi¬ 
cation  of  UD)  then  lead  to  the  different  QDPT  formalisms 
discussed  in  the  next  three  sections. 

Different  choices  of  V  produce  different  model  Ham¬ 
iltonians,  and,  while  they  give  the  same  infinite  order 
eigenvalues  and  eigenfunctions  [u a,  Eq.  (14)),  their 
truncated  (finite-order)  results  are  not  generally  equiv¬ 
alent.  Thus,  the  choice  of  subsidiary  conditions  to  com¬ 
plete  the  specification  of  U  may  affect  the  rate  of  con¬ 
vergence. 

It  is  well  known,  of  course,  that  once  UD  has  been 
specified,  Eq.  (22)  does  provide  an  explicit  equation 
for  Ux'  through  the  resolvent  operator 


Q  V  I  »K»I 

-H0'V«a-e( 


I  Q)  =  -Ra’  [#o>  Ux  11 «' 
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III.  CANONICAL  VAN  VLECK  PERTURBATION 
THEORY 

The  unitarity  of  the  decoupling  operator  U  In  VVPT  can 
be  ensured  by  expressing  it  in  exponential  form35,4* 

U=eG  ,  (26) 

with  O'  an  anti -Her mitian  operator 

G  =  -G’  .  (27) 

The  canonical  form  of  VVPT  (compare  Klein, 5  Jdrgen- 
sen, 42  and  Brandow4)  is  obtained  by  completing  the  spe¬ 
cification  of  U  through 

G=GX,  Gd=0.  (28) 

The  decoupling  operator  and  the  transformed  Hamilto¬ 
nian  obtained  from  this  condition  will  be  denoted  by  Uc 
and  X’c  =  F0  +  Wc ,  respectively. 

In  order  to  obtain  a  compact  formalism  we  shall  use  a 
superoperator  notation,  38,38  in  which  with  any  operator 
•A  we  associate  a  superoperator  A  defined  by50 

AX  =  [X,A\  =  XA  -AX  (29) 

(where  X  is  any  operator).  Positive  powers  of  A  pro¬ 
duce  repeated  commutators 

a2x  =  [[x,a|,  Aj  ,  (30) 

etc.,  and  the  zero  power  is  the  identity  superoperator 

A°.V  =  X  .  (31) 

As  noted  by  Primas,  38,38  this  allows  a  compact  repre¬ 
sentation  of  the  Baker-Campbell-Hausdorff  expansion 


Bernoulli  numbers  B (see,  for  example,  Abramowitz 
and  Stegun51): 


(39) 


(Other  notations  for  the  Bernoulli  numbers  are  also  in 
use;  see,  for  example,  Jolley.52)  The  first  few  coeffi¬ 
cients  in  the  series  are 


r®  - 1,  c,  =  3,  c*  =  - Cj- gf ft  C4  -  —  ihs*  *  *  *  • 

(40) 

Equation  (38)  can  be  expanded  order  by  order: 
[H0,GU>]  =  -VX  , 

[tf0,G<2l]=-(VD,G"’l  , 

[tf0,  G(3J  j  =  -  [Vfl,  G«M  -  Cn  G'1’]  , 

[H0,G(4>l  =  -[V0,Gl3']-i{t[V„Gl3>],G‘2'] 

+  [[Y*,G<2’l,Ga,J}  , 

[^o>  Gi5>  i  =  -  [V0,  Gl4,j  -  j{[[Vx,  Gll> ),  G,3>] 

+  [[Yjr,  Gl2)],  G<21]  +  [[ Vjf ,  G<31],  G(l>]} 

^[[[[Yx,G«'],G<IM,G<1M,G(3'].  (4D 

etc.,  and  converted  into  explicit  equations  for  the  G("’ 
using  the  resolvent  formalism  [Eqs.  (24)  and  (25)]. 
(Note  that  G,0’  =  0.) 


The  transformed  Hamiltonian  Xc  =  (x^)D  can  also  be 
obtained  compactly  in  terms  of  the  hyperbolic  functions 
of  the  superoperator 

=  coshG  Hd  +  sinhG  Hx 


e-cHec  =H+[H,G}+-~  [[//,  G),  G]  +  •  •  * 

=  =  ■  (32) 

Zo 

It  is  convenient  to  partition  this  expansion  into  even  and 
odd  functions  of  G: 

eG  =  coshG* sinhG  .  (33) 

Then,  if  G  satisfies  Eq.  (28),  we  find  that 

(X’C)D  =  (eG  H)d  =  coshG  H0  -  sinhG  tix  ,  (34) 

(X’c  )x  =  (eG  H)x  =  cosh(5  Hx  *  s  inhG  Hn  .  (35) 

The  decoupling  condition  (12)  can  now  be  written  in  the 
form 

0  =  PCc)x  =  coshG  Hx  *  [hd,  G 1  ,  (36) 

(note  that  the  superoperat  function  G'1  sinhG  involves 
no  inverse  powers  of  G  in  its  power  series  expansion). 
We  thus  obtain  a  commutation  relation  for  G: 

[Md,G\-~G cothG Hx  ,  (37) 

or,  noting  Eqs.  (9)  and  (10), 

[Ho,  G]  =  -  [V’D,  G  |  -  G  cothG  Vx 

m 

=  -{Vd,G\ -H^G^Vx,  (38) 

*r*0 

where  the  power  series  coefficients  r„  are  related  to  the 


=  Hd  +  [flD,  c  ]  +  sinhG  Hx 

=  tiD  -  (coshG  -  1)  cothG  Hx  -  sinhG  Hx 
=HD-  cschG(cosh2G  -  coshG  -  sinh2G)  Hx 
='HD  -cschG(l  -cosh G)HX 

=  HD  +  tanh(|G)Hx  ,  (42) 

Wc  =  VD  +  tanh(j  G)  Vx 


G)*‘XVX  , 


(43) 


where  the  power  series  coefficients  /„  are  also  related 
to  the  Bernoulli  numbers 


2>’»«(2»’»2.  i) 

,=  (2n  +  2)!  2 


(44) 


The  first  few  coefficients  in  this  series  are 

/0  =  1,  I1-  —  3,  fg-j },  h  =  ^ —  5§§T»  •••  •  (45) 

The  order-by-order  computation  of  H'c  (and  of  Vc) 
easily  follows: 

H'c"  =  1'd  • 
tf<2,=l[Vx.G",|, 

Hr=Urv.G<3M-WVG,,,|.G(,'].G,,'|. 
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-  [[[v„  G«M.  G«M,  G<,>]+[[[V’X,GW>],  G,l>],  G<uj}  , 

(46) 

etc.  We  observe  that  Wc*  depends  on  Gl"‘,  m  =  1,2, , 
n  -  1.  The  operator  Uc,  and  thus  the  perturbed  model 
functions  (13),  ca'  be  obtained  from  Eq.  (26)  to  the  same 
order  as  G. 

ft  is  clear  that  the  canonical  VVPT  formalism  is  en¬ 
tirely  expressible  within  the  domain  of  a  Lie  algebra, 
and  should  thus  lead  to  fully  linked  diagrammatic  expan¬ 
sions  which  satisfy  the  extensivity  requirement  order  by 
order. 


IV.  THE  KIRTMAN-CERTAIN-HIRSCHFELDER  FORM 
OF  VVPT 

The  canonical  form  of  WPT  was  derived  from  the 
choice  Gd=  0.  In  some  earlier  work  on  the  van  Vleck 
approach,53  and  more  recently  in  the  work  of  Kirtman43'44 
and  of  Certain  and  Hirschf  elder45, 48  (KCH),  a  different 
choice  was  made  to  complete  the  specification  of  the  de¬ 
coupling  operator.  It  does  not  appear  possible  to  spe¬ 
cify  this  choice  in  terms  of  a  single  order-independent 
equation,  so  it  is  stated  as  a  condition  connecting  dif¬ 
ferent  orders  of  V. 

The  unitarity  condition  of  U  through  any  order  n>  0 
can  be  stated  in  the  form 

n 

£  C"",,C/ta'",  =  0  (n>0)  .  (47) 

m«0 

The  sum  in  Eq.  (47)  can  be  split  into  two  parts 

+  V1  (1-4  i2*.„)  Ulm>  ,  (48) 

m»0 

where  [xj  denotes  the  integral  part  of  x.  The  KCH  ortho¬ 
normalization  condition  for  the  perturbed  model  func¬ 
tions  (13)  can  be  obtained  by  requiring  that  the  diagonal 
blocks  vanish  for  each  of  the  two  sums  on  the  rhs  of 
Eq.  (48)  separately.  Denoting  the  resulting  decoupling 
operator  and  transformed  Hamiltonian  by  UK  and  XK 
=  H0~  WK,  respectively,  we  get48 

<£•'*%  =  -  <1-5  5 <n>0).  (49) 

Together  with  the  uecoupling  condition  (12),  or  with  Eqs. 
(22)  and  (23),  this  completely  specifies  UK.  It  is  pos¬ 
sible  to  write  UK  in  exponential  form 

VK  =eK  ,  (50) 

with  K  =  -K\  but  there  appears  to  be  no  advantage  for 
this  form  in  the  KCH  treatment.43-48 


By  repeated  use  of  this  relation,  together  with  Eq. 
(49),  it  is  possible  to  convert  Eq.  (23)  into  a  form  in 
which  and  are  expressed  in  terms  of  U%*‘, 

m  =  1,  2, . . . ,  n.  The  resulting  equations  can  be  written 
in  the  form55 


W}?"  ={  {Up’1”  VU?>  +  u£"  VUX"U)D 


I»/2) 


(1  -5  60>.  -  4  6»,J  {W‘"‘(C/‘'— Uf-% 


+  (l/?-*'t{/£-*,*‘,)D  Wf'}  , 
WiK*“l>  =  (Up ,f  VU^) o 


(52) 


it  Un*l)A) 
ij  «0 


+  (Ujr*"ut,-m'l,",)D  Wf'}  ,  (53) 

and  represent  an  extension  of  Wigner’s  (2 n  +  1)  rule  to 
quasidegenerate  perturbation  theory.  This  should  result 
in  computational  economies,  and  is  the  principal  advan¬ 
tage  of  the  KCH  formalism.  43-48  It  does  not  appear  pos¬ 
sible  to  obtain  a  similar  result  for  the  other  QDPT  for¬ 
malisms. 


The  unitary  transformation  which  connects  UK  and  Vc 
can  best  be  expressed  in  terms  of  an  operator  .17,  de¬ 
fined  by 

UK -e“  =eae>l  =UC  e“  .  (54) 

This  operator  is  block  diagonal  {.Mx  =0)  and  anti-Her- 
mitian,  and  tedious  order-by-order  comparison  of  the 
equations  for  UK  and  Uc  shows  that 

,W<01  =M(1>  =Ma>  =0  , 

.Wl3>  =4  [Gll),  G,2>]  , 

,W(4>  =4  [G<u,  G131)  , 

m<5,=4  [g<1>,  g<4>] +4  [g<2,,g<3>] 

+  4  [G<*\(G<1,)3]-4G,l,[G<1\G(3,]G11’  ,  (55) 

etc.  No  simple,  order-independent  equation  for  .17  could 
be  found. 

A  recent  diagrammatic  treatment44  of  the  KCH  formal¬ 
ism  finds  that  the  expressions  for  Wjp'  can  be  written  in 
fully  linked  form  through  fourth  order,  but  that  some 
unlinked  diagrams  remain  in  fifth  order.  It  appears 
likely  that  this  is  related  to  the  lack  of  an  order-inde¬ 
pendent  equation  for  l'K  (or,  equivalently,  for  K  or  .17), 
since  it  can  be  shown  that  any  formalism  in  which  V  can 
be  expressed  in  the  exponential  form  (32)  satisfies  the 
separability  theorem28’8  and  is  thus  fully  linked,  order 
by  order,  provided  general  order-independent  equations 
for  the  exponential  operator  exist. 


V.  INTERMEDIATELY  NORMALIZED  QDPT 


Using  Eqs.  (16)  and  (20),  we  can  derive  the  recursion 
relation54 

«♦! 

T  yu<n>  1)  +  ^  £r(n> 

*«0 

(51) 

*•0 


The  most  widely  used  form  of  QDPT  is  based  on  in¬ 
termediate  normalization  and  leads  to  a  non-Hermitian 
model  Hamiltonian.2-10  We  shall  denote  the  decoupling 
operator  (wave  operator)  and  the  transformed  Hamilto¬ 
nian  obtained  from  the  intermediate  normalization  con¬ 
dition  by  U,  and  JC,  =H0*  \Vt,  respectively.  This  condi¬ 
tion  takes  the  simple  form 
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(V,)0  =  l  .  (56) 

It  is  convenient  to  define  an  operator5* 

X  =  (U,)x  =  a,  -  1  .  (57) 

Substituting  Eqs.  (56)  and  (57)  into  Eq.  (19)  and  then 
Eq.  (18)  provides  a  very  simple  direct  derivation  of  the 
equations  determining  X  and  W,: 

W,*VD  +  VxX,  (58) 

[H0,X\  =  ~  Vx  +  XVxX  .  (59) 

We  also  find  from  Eq.  (59)  that  X  is  anti-Hermitian: 

X  =  -X'  .  (60) 


Order -by-order  equations  for  X  and  follow  directly 
from  Eqs.  (58)  and  (59): 

n-2 

[tf„,  Y<n,]  =  -  [l,0,.Vto'u|+  21  Ar<”"  vx  A'0""'11 

m*  1 
rt^  2 

=  +  (”>2),  (61) 

=  v0  , 

w}n,  =  vx  A0"1’  (n  >  1)  .  (62) 

Comparing  the  recursion  relations  (61)  for  X  with 
those  for  G  [Eq.  (41)  j,  we  find  that 

G<1>  =Xlx>,  Gil>  =X<Z1  .  (63) 

A  more  tedious  comparison  shows  that 

G,3,=Xa>*j(X<l>)3  , 

G'*’  =X14’  +  U(X(,,)iX<i> +  X<1>X<z,X,v  +x<i>(x°>)2}  ■ 

(64) 

In  fact,  these  are  special  cases  of  the  general  relation¬ 
ship  between  G  and  X,  proved  in  the  Appendix, 

*  . 

G  =  arctanhY  =  - - r X 2"“  ,  (65) 

or 

X  -  tanhG  =  /„  G**1  (66) 

n*0 

[see  r.’S.  (44)  and  (45)  |.  Since  arctanlur  =  s  in[(l  *  .v)/ 

( 1  -  v)  1,  we  have57 


=  (1  -,Y)(1  -XlY'lz 

=  Vl(l’]VsYin  .  (67) 

As  seen  in  Sec.  II,  V)  U,  is  the  overlap  matrix  for  the 
perturbed  functions  l',\  I).  Thus,  as  noted  by  Brandow4 
(see  also  des  Cloizeaux, 37  Klein, 5  Kvasnicka,3  and 
Levy10),  the  canonical  van  Vleck  functions  Uc  I  a)  may 
be  viewed  as  the  result  of  the  symmetric  orthonormal- 
lzation5*  of  the  U,  i  a }  functions.  The  connection  between 
the  corresponding  model  Hamiltonians  is  then 


TCc  =ec’ff  ea 

=  (U'lUIYlnU)HVl(U]U,rin 

=  (f J]U,YlnUj  U,  X,(V)  V,  Y'n 

=  W,,U1)l/iK,(u;  a,)-''* 

=  (1-X*)U3X,(1  -XzYui  ,  (68) 

or57 

30/  =  coshG  Xc  sechG  .  (69) 

An  equation  for  the  direct  determination  of  the  Hermitian 
model  Hamiltonian  JCc  in  terms  of  X  instead  of  G  is  given 
in  the  Appendix. 

VI.  DISCUSSION 

Three  forms  of  quasidegenerate  perturbation  theory 
have  been  examined  here  in  terms  of  a  common  general 
formulation  based  on  a  similarity  transformation  of  the 
Hamiltonian.  This  approach,  together  with  the  sym¬ 
metric  treatment  of  the  P  and  Q  subspaces  represented 
by  the  D-X  (block-diagonal/off-diagonal)  partitioning  of 
operators,  has  enabled  simple  and  direct  derivations  of 
the  relevant  equations  and  has  clearly  brought  out  the 
relationships  between  the  different  QDPT  forms.  This 
derivation  bypasses  the  complications  of  the  usual  de¬ 
velopment  of  the  intermediate  normalization  formalism 
which  often  involves  iterative  removal  of  energy  depen¬ 
dence  from  the  denominators  and  treatment  of  zero- 
order  energy  differences  as  perturbations  (see,  for  ex¬ 
ample,  Kvasnicka, 8  but  see  also  his  more  direct  “al¬ 
gebraic  theory”  in  Sec.  Ill C>.  The  use  of  the  D-X  par¬ 
titioning  also  allows  us  to  avoid  the  extensive  use  of  pro¬ 
jection  operators,  which  tends  to  obscure  the  formalism 
(see,  for  example,  Klein5). 

Except  in  the  Kirtman-Certain-Hirschfelder  formal¬ 
ism,  43-4*  in  which  order-independent  formulation  is  not 
possible,  the  equations  for  the  wave  (or  decoupling)  op¬ 
erator  and  for  the  transformed  Hamiltonian  have  been 
obtained  initially  in  order -independent  (implicit)  form, 
from  which  recursive  order -by -order  equations  for  all 
orders  easily  follow.  The  use  of  a  superoperator  nota¬ 
tion  and  hyperbolic  functions  of  the  superoperator  has 
allowed  a  very  compact  and  simple  derivation  of  the 
general  equations  for  the  canonical  VVPT  formalism. 

The  infinite-order  model  Hamiltonians  of  the  different 
formalisms  are  related  by  similarity  transformations, 
but  when  these  model  Hamiltonians  are  truncated  at  a 
finite  order,  the  similarity  transformation  relationships 
are  no  longer  exact.  Thus,  they  would  generally  pro¬ 
duce  different  results  order  by  order,  with  possibly  dif¬ 
ferent  convergence  characteristics.  The  selection  of 
the  method  to  be  used  would  thus  be  governed  by  con¬ 
vergence  behavior  as  well  as  by  computational  consid¬ 
erations.4  In  this  regard  the  KCH  formalism  has  the 
advantage  of  the  (2«-  1)  rule,50  while  the  other  forms 
benefit  from  having  fully  linked  diagrammatic  expan¬ 
sions  in  all  orders.  The  intermediate  normalization 
form  (with  or  without  Hermitization)  appears  to  be  the 
most  convenient  for  applications  involving  infinite  order 
partial  summations.4 

Equation  (A5)  of  the  Appendix  provides  a  hybrid  ap- 
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proach  for  QDPT  calculations,  using  the  intermediate 
normalization  formalism  for  the  determination  of  the 
wave  operator  U,  =  1 +  X,  but  obtaining  the  Hermitian 
model  Hamiltonian  =H0+WC.  This  type  of  approach 
has  been  advocated  by  Brandow. * 
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APPENDIX 

In  this  appendix  we  derive  Eq.  (66),  connecting  the 
canonical  van  Vleck  and  the  intermediately  normalized 
forms  of  the  wave  operator.  We  also  present  an  equa¬ 
tion  expressing  jCc -Hn  (or,  equivalently,  Wc  -  VD)  in 
terms  of  Vx  and  X  directly. 

From  Eqs.  (11)  and  (26)  we  have 

H  =eaXc  e'a 

=  (coshG  +  sinhG)  JQ;(coshG  -  sinhG)  .  (Al) 

Taking  the  block -diagonal  and  block  off-diagonal  parts 
of  this  equation  separately,  we  find 

Hd  =  coshG 3Cc  coshG-  stnhG  sinhG  ,  (A2) 

Vx  =  -  coshG  sinhG  +  sinhG  5C<;  coshG  .  (A3) 

Therefore, 

[f/0,  tanhG]  =  (coshG  jCc  coshG  -  sinhG  3cc  sinhG)  tanhG 

-  tanhG(coshG3Cc  coshG  -  sinhG jQ;  sinhG) 

=  coshGJCc  sinhG  -  sinhG  JCc  coshG 

-  tanhG(coshGJCc  sinhG  -  sinhG  a’c  coshG)  tanhG 

=  -  Vx  -  tanhG  Vx  tanhG  .  (A4) 

Thus,  tanhG  satisfies  the  same  commutation  equation  as 
X  [Eq.  (59)  j.  Since  also  (tanhG)o=0  [this  follows  from 
Eq.  (28)  and  the  fact  that  tanhG  is  an  odd  function  of  G], 
tanhG  must  be  identical  to  the  X  operator. 


The  equation  referred  to  above  for  \VC  -  VD  is  obtained 
by  a  lengthy  and  tedious  sequence  of  algebraic  manipula¬ 
tions,  beginning  with  Eq.  (68)  and  using  Eqs.  (58)  and 
(59)  and  the  properties  of  the  binomial  expansions  of 
(1  -  x)*1/s.  The  result  can  be  expressed  in  any  of  the 
forms 


WC‘VD~i{Vx,X\ 


2  (n  *  2»i)(«  *  2ih  +  1) 


=  Vn 


c  {.Y2”[VX ,  .V2"’1 1 X *  A'2"*1  [vx ,  A'2"'1 1  A'2"*1 

dvx,x\ 


lyy _ n- - 

2  TTl  («  *  2 »«)(n  +  2m  ■ 


x[VxX  ~XVx,X*\Xln 

•Vo*HVx.X\ 


1  vv 


(n  ♦*!>()> 


X*"(VxX  +  XVx)r/nX*  , 
(A5) 


where  ' 


is  the  coefficient  of  V  in  the  binomial  power  series  for 
(1  -  .r)‘1/2.  The  first  two  terms  in  these  expressions 
represent  the  Hermitian  average 

(W/)M  =  i(W/  +  lV,,)  =  V0  +  i[Vx,Aj  .  (A7) 

The  additional  terms  contribute  to  Wc  beginning  in  fourth 
order.  The  third  order  form  of  this  model  Hamiltonian 
has  been  used  by  Freed  and  co-workerss0  in  their  work 
on  the  effective  valence  shell  Hamiltonian. 
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Multidimensional  many-body  theory:  Diagrammatic 
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A  size-extensive  multidimensional  many-body  theory  is  developed  from  an  order-expanded  van  Vleck 
transformation.  This  provides  an  effective  Hamiltonian  in  a  model  space  consisting  of  a  set  of  determinants 
whose  zeroth-order  energies  may  be  nondegenerate.  Expressions  for  the  effective  Hamiltonian  in  terms  of  the 
perturbation  and  a  set  of  resolvents  generalized  from  the  Rayleigh-Schrodinger  form  are  given.  Perturbative 
evaluation  of  the  resultant  formulas  via  diagrammatic  expansion  is  illustrated  and  discussed.  The  diagrams 
required  through  second  order  for  a  model  space  consisting  of  a  Hartree-Fock  solution  plus  selected  singly 
and  doubly  excited  determinants  are  presented,  and  their  relation  to  those  employed  in  the  method  of  Hose 
and  Kaldor  is  discussed. 


I.  INTRODUCTION 

Two  problems  of  current  theoretical  interest  are  (1) 
the  calculation  of  highly  correlated  excited  state  ener¬ 
gies  and  potential  surfaces  for  molecules  and  (2)  the 
description  of  ground  states  of  multiply-bonded  or 
open-shell  species.  The  first  ab  initio  method  used 
for  these  purposes  was  that  of  configuration  interac¬ 
tion  (Cl). 1  Original  applications  involved  only  a  few 
selected  configurations,  but  today  they  routinely  include 
all  single  and  double  excitations  relative  to  a  given 
configuration  or  set  of  reference  configurations. 

In  recent  years,  the  Cl  approach  has  been  joined  by 
many-body  perturbation  theories  (MBPT)2-5  [based  on 
Rayleigh-SchrOdinger  perturbation  theory  (RSPT),  dia¬ 
grammatic  expansions,  and  the  linked  cluster  theorem 
of  Goldstone3],  including  iminite-order  coupled-cluster 
techniques. 3-10  MBPT  calculations  based  on  a  single 
determinant  have  become  routine  for  ground  states  of 
many  systems,3’11  and  their  usefulness  for  a  wide  vari¬ 
ety  of  molecular  problems  has  been  demonstrated. 12-15 

In  MBPT,  it  has  proved  convenient  to  adopt  a  dia¬ 
grammatic  formulation  since,  just  as  it  is  not  feasible 
to  solve  the  full  Cl  problem  in  a  basis  of  sufficient  size 
to  quantitatively  study  interesting  molecular  systems, 
the  RSPT  energy  expansion  must  be  truncated.  The 
linked  diagram  expansion  provides  a  useful  subdivision 
of  the  terms  to  be  evaluated  or  neglected.  Even  so, 
complete  evaluation  of  fourth  order  involves  terms  from 
triple  and  quadruple  excitations  (relative  to  the  Har¬ 
tree-Fock  ground  state).  Various  fourth-order  ap¬ 
proximations  (or  particularly  the  analogous  Cl  treat¬ 
ments)  are  still  formidable  problems.  On  the  other 
hand,  stopping  with  third-order  RSPT  does  not  incluJe 
any  effects  of  single  excitations  or  higher-than-double 
excitations.  Furthermore,  any  less-than-full  Cl  treat¬ 
ment  or  nondiaerammatic  RSPT  approximation,  which 
is  not  complete  to  a  given  order,  will  suffer  from  size 
inextensivitv  (also  called  size  inconsistency  or  un¬ 
linked  cluster)  errors;  i.e. ,  the  energy  is  not  propor¬ 
tional  to  the  number  of  noninteracting  subsystems, 
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The  connection  of  RSPT  to  diagrammatic  perturbation 
theory  and  the  linked  cluster  theorem  of  Goldstone 
provides  a  solution  to  these  problems.  Each  diagram 
possesses  the  important  property  of  size  extensivity. 
This  means  that  approximations  consisting  of  evaluation 
of  subsets  of  the  full  set  of  diagrams  for  a  given  order 
will  also  have  this  property  as  long  as  complete  dia¬ 
grams  are  included.  Furthermore,  diagrammatic 
studies  have  enabled  the  observation  that  often  the  ma¬ 
jor  function  (in  fourth  order)  of  including  the  Cl  quadru¬ 
ple  excitations  is  merely  to  cancel  the  spurious  (un¬ 
linked)  parts  of  the  effect  of  double  excitations  in  Cl, 
which  are  not  included  in  the  linked  cluster  expan¬ 
sion.  10,17  Most  importantly,  diagrammatic  perturbation 
theory  has  led  to  formulations  for  the  further  effects 
of  quadruple  excitations  which  are  feasible  to  calcu¬ 
late.  10’11 

Although  one-dimensional  MBPT  has  been  shown  to 
provide  highly  accurate  binding  and  activation  ener¬ 
gies, 18-20  potential  surfaces,21-25  and  spectroscopic 
parameters,  such  as  the  quartic  force  field  of  the  water 
molecule,13  it  is  not  universally  applicable,  particularly 
when  exact  or  near  degeneracies  are  encountered  in  a 
molecular  problem.  For  example,  even  with  unre¬ 
stricted  Hartree-Fock  reference  functions,  it  is  usual¬ 
ly  insufficient  for  potential  surfaces  involving  the  dis¬ 
sociation  of  multiple  bonds, 11  and  it  is  not  ordinarily 
applicable  to  excited  states  other  than  the  lowest  state 
of  each  symmetry.  The  study  of  these  problems  neces¬ 
sitates  a  multidimensional  treatment. 

Several  multidimensional  perturbative  methods  have 
been  proposed.  The  completely  degenerate  case  was 
considered  by  Bloch  and  HoroWitz, 2i  and  later  by 
Morita, 2:  The  first  thorough  quasidegenerate  linked- 
diagram  approach  was  formulated  by  Brandow.28 
Lindgren  presented  a  direct  and  particularly  trans¬ 
parent  development  of  the  theory,  which  he  then  ex¬ 
tended  to  a  multidimensional  generalization  of  the 
infinite-order  coupled-cluster  approach.29  More  recent¬ 
ly,  Hose  and  Kaldor  (HK),30  Mukherjee  el  a!.  ,31  and 
others32-35  have  contributed  to  this  field  of  study.  The 
results  of  Freed  and  co-workers, 35  Kaldor,33  Hegarty 
and  Robb,37  and  HK38  indicate  the  promise  of  these 
theories. 
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TABLE  l.  Characterization  of  methods. 


Method 

a 

b 

c 

d 

B  randow* 

i 

N 

c 

S(C) 

F  reed‘ 

i 

S 

c 

N 

Hose  and  Kaldor*’11 

i 

N 

G 

T 

Llndgren1 

D 

N 

C 

S(C> 

Present 

D 

H 

G 

S(HF) 

‘Direct  (D)  or  indirect  (I)  formalism  for  denominators. 
“Hermitean  (H).  non-Hermitean  (N),  or  symmetrized  (S)  effec¬ 
tive  Hamiltonian  formalism. 

‘General  (G)  or  complete  (C)  model  space  implementation. 
Single  reference  state  (S).  such  as  core  state  (C)  or  Hartree- 
Fock  determinant  (HF);  or  transition  (I)  diagrammatic  nota¬ 
tion;  or  nondiagrammatic  implementation  (N), 

'Reference  2d. 

‘Reference  35. 

“Reference  30. 

"Reference  38. 

‘Reference  29. 


The  relationships  among  these  methods  can  be  ex¬ 
amined  on  both  formalistic  and  implementational  levels. 
One  formal  difference  involves  whether  the  energy 
denominators  are  obtained  directly  (D)  or  indirectly 
(I,  employing  a  degeneracy-breaking  perturbation  and 
resumming).  Also,  formalisms  can  lead  to  Hermitean 
(H)  or  non-Hermitean  (N)  effective  Hamiltonians.  Some 
of  the  former  are  explicitly  Hermitean,  while  others 
are  symmetrized  by  brute  force  (S).  This  aspect  is 
analogous  to  the  choice  of  fully  or  intermediately  nor¬ 
malized  perturbed  model  functions.  It  has  been 
shown39,40  that  the  canonical  van  Vleck  functions  can 
be  viewed  as  the  symmetrically  orthogonalized  version 
of  the  intermediately  normalized  functions. 

On  the  implementational  level,  the  basic  difference 
is  whether  a  complete  valence  model  space  is  required 
(C)  or  the  model  space  is  allowed  to  be  general  (G, 
or  incomplete).  A  related  matter  for  the  specification 
of  diagrams  is  the  choice  of  reference  state(s)  from 
which  particles  and  holes  are  defined.  The  options  here 
are  the  use  of  a  single  state  (S),  which  could  be  a  core 
state  with  less  than  N  orbitals  or  a  given  state  (e.  g.  , 
the  Hartree-Fock  determinant,  as  in  the  present  meth¬ 
od),  or  the  use  of  several  different  states  [as  in  the  HK 
transition  (T)  method].  Each  of  the  various  methods  of¬ 
fers  a  combination  of  these  features.  Table  I  sum¬ 
marizes  the  characteristics  of  several  approaches,  in¬ 
cluding  the  present  one,  using  the  notations  introduced 
above. 

Brandon  's-’3  non-Hermitean  wave  operator  approach 
uses  a  core  reference  and  a  complete  model  space.  It 
expands  the  Brillouin-Wigner  (BW)  denominators,  rela¬ 
tive  to  a  reference  energy,  and  incorporates  a  de¬ 
generacy-breaking  perturbation,  from  which  a  Ruv- 
leigh-SchrSdinger  size-extensive  expansion  is  re¬ 
gained  via  resummation  techniques.  Although  it  is 
equivalent  to  Brandow's  approach,  Lindgren’s  for¬ 
malism3’'  uses  a  generalization  of  the  degenerate  Block- 
Horowitz3  method  to  dispense  with  the  energy  depen¬ 


dence  of  the  effective  Hamiltonian  initially,  and  yield 
a  non-Hermitean  order -by -order  perturbation  expan¬ 
sion  for  the  wave  operator,  without  considering  BW 
denominators.  In  addition,  when  the  cluster  operator 
es  is  put  into  normal-product  form,  a  convenient  de¬ 
velopment  of  multidimensional  coupled-cluster  theory 
emerges. 

The  approach  of  Hose  and  Kaldor30  differs  in  two 
significant  ways  from  the  earlier  methods.  First,  it 
permits  a  general  model  space  rather  than  a  complete 
space.  This  flexibility  is  important  for  reducing  com¬ 
putation  time  and  for  handling  the  problem  of  intruder 
states.  Also,  it  specifies  each  matrix  element  in  per¬ 
turbation  theory  with  respect  to  its  ket  vector  (e.g. ,  the 
ket  vector  of  Hai  is  13)),  rather  than  to  a  single  core 
state.  This  approach  introduces  some  cancellations  that 
are  not  present  in  the  Brandow  development.  These 
authors  handle  their  energy  denominators  as  Brandow 
does,  in  a  quasidegenerate  fashion. 

The  approach  proposed  in  this  paper  appears  to  have 
advantages  over  these  earlier  techniques.  We  believe 
it  is  necessary  that  multidimensional  methods  be  size- 
extensive,  permit  arbitrary  choices  of  model  func¬ 
tions,  and  still  retain  clearly  defined  hole  and  particle 
labels  so  that  categories  of  higher  excitations  may  be 
easily  identified.  A  method  that  offers  these  elements 
and  is  manifestly  Hermitean  is  the  canonical  van  Vleck 
transformation.  40-42  Brandow  has  discussed  the  simi¬ 
larity  of  this  method  to  his  apprach, 39  but  he  has  not 
considered  the  generality  of  the  method  for  arbitrary 
model  spaces. 

Another  diagrammatic  (noncanonical)  van  Vleck  ap¬ 
proach  has  been  proposed  by  Kirtman.33  This  formalism 
(referred  to  as  KCH  in  Ref.  40)  has  the  advantage  of  a 
2«-  1  rule;  i.e. ,  The  (2»-  l)th-order  effective  Hamil¬ 
tonian  can  be  expressed  in  terms  of  the  nth-order  per¬ 
turbed  functions.  His  results  are  size-extensive  through 
fourth  order,  but  not  in  higher  orders,  in  which  onlv  a 
‘‘quasicancellation”  of  unlinked  diagrams  occurs  unless 
the  model  states  are  exactly  degenerate. 

In  the  following,  we  implement  the  canonical  van  Vleck 
transformation  formalism  for  nondegenerate  MBPT  in 
an  arbitrary  multidimensional  model  space.  We  present 
an  example  of  the  diagrammatic  expansion  through  sec¬ 
ond  order.  The  advantages  of  this  method  are:  ( 1 ) The 
method  yields  a  Hermitean  effective  Hamiltonian  and. 
thus,  a  set  orthogonal  eigenvectors;  (2)  the  resultant 
energies  will  be  specified  bv  an  expansion  of  size- 
extensive’  diagrams;  (3)  in  our  diagrammatic  expansion, 
we  observe  cancellations  that  do  not  occur  in  the  other 
methods;  14)  because  we  construct  all  diagrams  wun 
reference  to  a  Hartree-Fock  determinant  UKF  for  ar. 
open-shell  problem),  higher-excitation  cont.-ibjtkv'.- 
are  easily  identified,  and  useful  approximations  for 
actual  compulations  mav  be  more  readilv  apparent: 

(5)  a  single  evaluation  and  diagonaliza'icn  of  the  ef¬ 
fective  Hamiltonian  will  vield  a  sei  oi  neariv  equiva- 
lentlv  treated  slates  (bias  resulting  :reir.  the  choice  ot 
orbitals,  a  common  problem  shared  bv  all  mullisiatt 
methods,  is  not  considered  here):  and  (6'  implementa¬ 
tion  within  a  general  model  space  of  selected  deter¬ 
minants  increases  the  practicality  of  applications. 
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II.  THE  CANONICAL  VAN  VLECK  (CWI  along  with  Eq.  (7),  and  analogous  terms  for  the  blocks 

TRANSFORMATION  of  H ,  V  and  U,  we  have 


Multidimensional  perturbation  theories  begin  from  the 
specification  of  a  model  space  spanned  by  functions 
{#,,  ysl,  d\.  If  the  *’s  are  orthonormal,  the  projectors 
P  for  the  model  space  and  Q  for  its  orthogonal  comple¬ 
ment  are  given  by 

i 

i>  =  £  !+rK*rl*r>',(*rl=£  |v><>  I  ■  (D 

r*l  r 

<?  =  £  !♦{> |  =  £  | i>< * |  .  (2) 

l>  i  * 

P  and  Q  are  idempotent  and  self-adjoint,  and  form  a 
resolution  of  the  identity.  The  objective  is  to  solve  for 
a  subset  of  eigenvalues  (and  eigenvectors)  of  the  Schro- 
dinger  equation 

HU  =  *E  (3) 

by  means  of  an  effective  Hamiltonian  equation  in  a  space 
of  reduced  dimension,  that  of  the  model  space 

3C.,fD  =  DE  .  (4) 

X*t  is  composed  of  matrix  elements  within  the  model 
space  of  a  transformed  operator 

X=UmtHU  ,  (5) 

D  is  a  matrix  of  eigenvectors,  and  E  is  a  diagonal  ma¬ 
trix  of  eigenvalues.  The  transformation  U  is  required 
to  be  unitary  and  to  decouple  the  P  and  Q  spaces43 

Uml  =  U*  ,  (6) 

<CX  sPXQ+QKP=Q  .  (7) 

The  exact  eigenfunctions,  if  needed,  are  given  by 

*  =  UP*D  .  (8) 

The  unitarity  of  U  results  in  full  normalization  since 

('H*)  =  Df(*|PUtUP*)D  =  l  .  (9) 

We  identify  P*D  as  the  perturbed  model  functions 

*°=P*D  =  t'f*.  (10) 

Using  Eqs.  (10)  and  (5)  in  the  form 

Hi'  =  UK  ,  U) 

Eq.  (4)  is  easily  derived  by  the  following  sequence  of 
manipulations  of  Eq.  (3): 

H*  =  * E  , 

HL'V’*=  l/l|,*E  , 

HU*°  =  C'*°E  , 

U  iC*J  =  6’*JE,  (12) 

=  *JE  , 

PXP* D  =  P$DE  , 

•P'PXP* D  =  DE  . 

Perturbative  expressions  for  U  and  may  be  obtained 
from  Eq.  (11).  Using  the  notations  of  Ref.  40 

<D  =p<p~Q<q,  (13) 

and  the  usual  separation  of  the  Hamiltonian  H  =  H{0)  *  V, 


Hz  —  Vx  ,  (14) 

Hd  —Hw+Vd  ,  (15) 

and  we  define'  W  by 

(16) 

Partitioning  Eq.  (11)  into  diagonal  and  off-diagonal 
blocks  leads  to 

[hw,us)=-vxud-vdux  +  uxw  ,  U7) 

and 

UDW  =  [Ha>,  UD\+VXUX  +  VDUD  .  (18) 


When  U  and  <C  (or  W)  are  expanded  in  orders  of  the  per¬ 
turbation  with 

U<0,=  1  ,  (19) 

3Cio»_H(oi  (  (20) 

(»>0)  ,  (21) 

these  equations  lead  to  explicit  recursive  equations  for 

U1"’  and  W4"4  (for  n>0): 


"-1 

[Hw,  Uif)]=-VxUo”'" -VDUjin-"  +  £  , 

mm  t 

(22) 

W*"1  =[HIC>,  (/<”> J  +  Ux0’x‘”-1'+  VDUinmt>~  £  . 

m-\ 

(23) 

Unitarity  is  ensured  by  defining 

V  =  f  ,  (24) 

with  C  being  anti-Hermitean 

C  =  -Gt  .  (25) 

It  is  then  convenient  to  obtain  ordered  equations  in  terms 
of  G  rather  than  V.  We  note  that  V  has  not  yet  been 
completely  specified.  This  is  accomplished  by  the  Kem¬ 
ble  condition, 11  requiring  that 

Go  =  0  .  (26) 


This  ensures  that  G  has  the  least  possible  effect  upon  the 
model  functions33  and  that  the  perturbed  model  eigen¬ 
vectors  optimally  resemble  the  exact  ones. 39,u 


From  Eqs.  (24)  and  (5),  it  follows1"'39  that 

<=e'aHec  =  We°)L  .  (27) 

The  subscript  L  means  that  all  effective  Hamiltonian 
matrix  elements  are  represented  entirely  in  terms  of 
“Unked”  diagrams.  In  general,  linked  diagrams  must 
not  contain  separate  closed  parts.  However,  diagrams 
consisting  solely  of  noninteracting  open  lines  and  a 
single  closed  part  (which  appear  on  the  diagonal)  are 
considered  linked.  Disconnected  open  diagrams  are 
also  permitted.  The  meaning  of  “open"  is  given  in 
Sec.  HI. 

As  is  shown  elsewhere40  the  expansions  of  G  =  G,  and 
<0  are  given  in  commutator  form  by 


J.  Chem.  Phvv,  Vol.  76,  No.  4,  15  February  1982 


14 


L.  T.  Redmon  and  R.  J.  Bartlett:  Multidimensional  many-body  theory 


1941 


[tf43’, 

G44'] 

l=-v*, 

[w43’. 

G421] 

=  ~  [ 1  D  > 

G("l 

. 

[Hia\ 

G13’] 

\  =  ~  (  > 

g‘2,| 

,  G41’]  , 

(Hl0\ 

Gu’l 

=  ~  [  I’d  > 

g43‘] 

-mu'x.G4*’ 

l,c42'! 

,  g42’ 

],  G411]}, 

[tf40’, 

G<s,l 

=  -[  V 

G44’] 

-H[(lx,G41) 

1,  G43,l 

*[\vXi 

,  G42’ 

1  G42>l  +  [[Vx> 

g43>1,  g1 

+  (i/45)[[[[rx><;‘,,l,  G4t>),  c41’],  g(U] 

(2! 

and 


k.'40,= 

H40’  , 

K.’411  = 

IV411  = 

I’d 

, 

^2,= 

IF421  = 

it 

lx. 

G4" 

•C43’  = 

VI’4  3  *  = 

it 

I'x. 

G421 

H.’441  = 

if44)  = 

it 

I’x- 

G43’ 

- (l/24)[trvx,  G"’],  G“'l,  G"*,!  , 

<c*s>  =  v^sl,=  HKx,  g44’] 

-  (I/24){([[PZ,  g“,1,  G“>],  G42’] 

+  [[[  ,  G4t>],  G‘2’l,  G41’] 

+  [[[KXlG‘2>])Gtl,})G<"]}.  (29) 

The  Hermiticity  of  the  expansion  for3C,f,  is  apparent, 
and  no  explicit  symmetrizing  procedure  need  be  applied. 
Unlike  Kirtman’s  generalized  van  Vleck  approach, 33  no 
powers  of  GU)  appear  within  the  commutators  of  Eq. 

(29). 15  Thus,  the  results  of  the  canonical  van  Vleck 
model  will  be  size-extensive  since  >C  is  represented  in 
the  domain  of  a  Lie  algebra.  u  Furthermore,  this  prop¬ 
erty  is  not  dependent  on  using  a  complete  reference 
space.  Thus,  this  treatment  incorporates  a  flexibility 
crucial  to  solving  the  multidimensional  many -body  prob¬ 
lem. 

The  explicit  diagrammatic  construction  of  the  effec¬ 
tive  Hamiltonian  matrix  elements  follows  conveniently, 
as  in  the  one-dimensional  MBPT  case,  from  expres¬ 
sions  given  in  terms  of  the  perturbation  V  and,  in  this 
case,  a  set  of  generalized  Rayleigh-SchrOdinger  resol¬ 
vents.  These  resolvents  arise  in  iC  from  the  commuta¬ 
tor  with  H <0’  in  the  equation  for  G.  The  expressions  for 
G<nl  are  realized  by  substituting  Eq.  (28)  into 


G* 

where 


- Z  Gx]py-pr[it",  cx]*h  . 

(30) 


Pr=  IrXvl  ,  (31) 

and  R is  a  generalized  Rayleigh-SchrCdinger  type 
resolvent  involving  the  zeroth-order  energies  of  the 
complementary  states  and  a  particular  model  state: 

*r=  E  .  02) 


Gm=  z  . 

>■ 

G‘2’=  £l*HV'D,  G{l’]Py-Pr[VD,  G(1,l , 

r 

Gl3,=  Z  G‘3,l-Htl'*,  Ga'l,  G“'])Pr 

r 

-Pr([v0,  g12,]  +  M[[^x.  G“'l.  Gu,!)p<<>'t , 

Gu,=  £  tfi‘0’((VD,  G<3)  ]  +  y  ([[  Vx ,  G,wl,  G<2'1 

r 

[[  l'x  >  G<2>],  G(  *'  ]))  Pr  -  PT 

>  (■  •  •  same  terms  •  •  ■  )Pj,01}  .  (33) 

These  expressions  are  sufficient  to  obtain  the  energy 
through  fifth  order.  The  equations  are  quite  lengthy  in 
tl.-  ;r  Rayleigh-SchrBdinger  form,  so  we  shall  give  the 
results  only  through  third  order: 

=  {4>B  I  , 

*iV  =  5  P**/  . 

=  I  j<*.  1  VR^  ’  VR‘a°  ’  V<t>B  > 

-  Z  <+jr*,K*,|KR‘4,/l'.,,r*|/ 

r 

■^(^J.vr^vr^v^) 

-  Z  <4».  I  V»i0>«i0>K<t.y>  <*„ ;  V4vj.  .  (34) 

r  1 

Each  order  (t3)  includes  renormalization  terms  simi¬ 
lar  to  those  resulting  from  the  bracketing  theorem  of 
ordinary  RSPT,  except  for  the  specific  combinations  of 
the  various  model  state  resolvents.  Through  third  or¬ 
der,  the  CVV  effective  Hamiltonian  of  Eq.  (34)  is  exactly 
equivalent  to  a  Hermitized  combination  of  the  Hamil¬ 
tonian  of  Lindgren25  or  HK.30  It  already  differs  from 
Kirtman’s  Hamiltonian  in  third  order. 

III.  DIAGRAMMATIC  EXPANSION 

In  the  previous  section,  we  obtained  Rayleigh-Schrb- 
dinger  type  formulas  for  the  order -by -order  determina¬ 
tion  of  an  effective  Hamiltonian  matrix,  which  when  dia¬ 
gonalized  will  yield  a  set  of  ground  and  excited  states. 

At  that  point,  the  type  of  model  functions  and  the  dia¬ 
grammatic  expansion  were  left  unspecified.  The  main 
objectives  of  a  multidimensional  treatment  are  (1)  to 
obtain  sets  of  uniformly  described  excited  states  and 
(2)  to  properly  correlate  ground  state  potential  surfaces, 
which  require  more  than  a  single  reference  function  be¬ 
cause  of  spin  symmetry  requirements,  breaking  of 
multiple  bonds,  or  open-shell  products,  etc.  We  will 
consider,  as  an  example,  a  model  space  consisting  of  a 
Hartree-Fock  (SCF)  determinant  plus  selected  singly 
and  doubly  excited  determinants. 

We  begin  by  specifying  the  separation  of  the  Hamil¬ 
tonian.  In  second  quantization,  the  electronic  Hamil¬ 
tonian  is  given  by 


These  energies  will  be  specified  in  the  following  section. 
Thus, 


»=Z  vov  +  i  Z 

N  ’  Hrl 


\ft\\ rs)p,^sr  , 


(35) 
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where 

**=  /  £  ^]x.«)dT,  .  (36) 

>/  rs)  =  /  x;(l)\,*(2)rii(X-/>12)Xr(l)x,(2)dr1dr,  , 

(37) 

and  we  have  designated  the  creation  and  annihilation 
operators  by  their  subscripts  alone.  In  terms  of  quan¬ 
tities  in  normal  product  form, 15  ff  is  given  by 


ff=ff*  +  (*»lff*#> , 

(38) 

where 

ff*  =  ff£  +  t*  , 

(39) 

^  =  23  +  23  ;t/()j‘v[/>,<7] . 

(40) 

!■*=  7  23  vpr/  rs/.VfpV/sH  . 

4  Pora 

(41) 

The  use  of  the  V*  potential  significantly  reduces  the 
number  of  diagrams  when  a  SCF  orbital  basis  is  em¬ 
ployed.  i;  Alphabetically,  («,&,...)  will  be  reserved 
for  virtual  orbitals  of  the  Hartree-Fock  solution, 

(i,  j, . . .)  will  be  for  those  orbitals  originally  occu¬ 
pied,  and  (/>,  q,  . .  . )  will  be  undesignated.  Particle  or 
hole  designation  is  always  defined  with  respect  to  4>0. 


The  SCF  orbitals  satisfy 

>‘N  +  £  <  pt  j  I  qt >  =  6*  ,  (42) 

t 

so  that 

<4.0jff!4>,)  =  £?F  .  (43) 

Thus,  we  have  the  separation 

ff  =  ff°+FN,  (44) 

with 

ff°  =  ff,:<  +  E?F  .  (45) 

The  model  states  are  defined  in  terms  of  creation 
and  annihilation  operators  working  on  the  Hartree-Fock 
determinant.  A  singly  excited  state  will  be  denoted 

'<fr‘;,  =  N[atil  !+o)  •  (46) 

Similarly,  a  doubly  excited  model  determinant  is  given 
by 

,  4’?/ )  =  -V [ ar  br  ji  J  j  4v  •  (47) 

Note  the  equivalences 

N(dff]  =  atf  ,  (48) 

N  }  =  tflfji  =  a  itf  j  ,  (49) 


resulting  from  the  designation  of  occupied  and  virtual 
labels. 

Diagrammatically,  each  model  function  is  represented 
by  a  diamond  vertex  with  appropriately  labeled  particle 
and  hole  lines  entering  and  leaving.  Thus,  any  vector 
originating  or  terminating  at  a  diamond  vertex  will 
represent  a  specific  orbital  rather  than  an  index  of 
summation.  These  lines  are  the  so-called  "open" 


TABLE  II.  Zeroth-order  effective  Hamiltonian 
matrix  elements. 

HFI  HF: 

<f>0)  =£o F 
H  F I SE  =  U 
H  F  I  DE  =  0 
SE  I  SE: 

SE 1 DE  =0 
DE  I  DE: 


lines.  A  diagram  containing  open  lines  may  also  be 
described  as  open.  The  conv  ntion  of  not  redefining 
hole  and  particle  labels  results  in  allowing  open  lines 
on  either  side  of  the  diagrams. 

A.  Zeroth  and  first  orders 

The  zeroth-order  effective  Hamiltonian 

■<iV  =  <*.;»?*•  >  +  6«»4F  (50) 

is  diagonal.  Matrix  elements  between  the  non-,  singly-, 
and  doubly -excited  determinants  (HF,  SE,  and  DE, 
respectively)  are  evaluated  using  the  generalized  Wick’s 
theorem. The  results  relative  to  £“F  for  the  chosen 
model  space  are  given  in  Table  II.  We  see  that  ff "  acts 
as  an  excitation  operator. 

The  first-order  contributions 

<V  =  <+.l^"+»>  (5D 

are  evaluated  similarly.  For  the  Hartree-Fock  deter¬ 

minant 

( 4>ol  FN$o) =  0  >  (52) 

because  of  the  normal  product  form  of  t'N.  In  ac¬ 
cordance  with  Brillouin’s  Theorem,  we  obtain 

(*ol  FN4>? >  =  (4>ol  T  23  \pq<  rs) 

*  Hr* 

*N[/>,qTs>'lN[(j,i]<l>0l>  =  O  ,  (53) 

because  the  two  normal  products  cannot  be  fully  con¬ 
tracted.  The  diagrammatic  analog  is  that  there  is  no 
way  to  properly  connect  the  Hartree-Fock  bra  (no 
lines),  the  perturbation  interaction  (four  lines),  and  a 
single  excitation  ket  (two  lines). 

The  first-order  interaction  between  the  Hartree-Fock 
determinant  and  a  double  excitation  from  it  is  given  by 

<4>0 ]  VrN  )  =  <4>0  t  23  (/"?  rs,'N[pV*H 

VN  [n’i’yi]  <t>o- 

1  V'  , 

—  7  \P*i  ts  <5^,,  5,^ 

’  Ptr> 

■=\ij'  ab)  ,  (54) 
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TABLE  01.  First-order  effective  Hamiltonian  matrix 
elements. 

HF  I  HF  =  0 
HFI SE  =0 
HF  I  DE: 

(*,|V 

SE  I  SE: 

SE  I  DE: 

(*}IFH  *?,)  <aUJb)  I1crf» -{phr} 

OE 1 DE: 

( *(J  I  F1*  )  =  a,*,,*,  6ljt)iu)  < (yXdc)  II  (abm) ) 

+  {6jyU  {ab  llcd)}  +  {phr} 


where  the  6'  notation  is  defined  by 

i  s  6>(  6„  ~  6>y  •  (55) 

Proceeding  as  in  these  examples,  we  find  the  results 
listed  in  Table  III,  where  the  notation 

6,1  „ » (X  m  I  |X>0  =  6,»<Xf  || XX)  -  6„  Qik I  |XX)  (56) 

has  been  used.  Similarly, 

6 (X (dc)\\(ab)X)  =  6tc (Xd 1 1 oX>  -  6 M(Xc 1 1 aX) 

-  5«  CXd  ||  &X>  +  e^CXc  1 1 6X)  . 

(57) 

Thus,  the  first  term  of  the  last  matrix  element  in  Table 
HI  represents  16  possible  terms,  all  of  which  have  the 
same  Hugenholtz  diagrammatic  skeleton. 

The  first-order  effective  Hamiltonian  energy  diagram 
forms  are  summarized  in  Fig.  1.  Note  that  terms  in¬ 
volved  in  the  parenthetic  notations  always  share  a  com¬ 
mon  form.  Forms  related  by  particle-hole  reversal, 
which  will  be  omitted  for  brevity  in  subsequent  figures, 
have  been  shown  explicitly. 

The  algebraic  expressions  can,  of  course,  be  derived 


from  the  diagrams.  As  usual,  it  is  necessary  to  ex¬ 
pand  the  (antisymmetrized)  Hugenholtz  diagrams  to  (anti- 
symmetrized)  Goldstone  form  in  order  to  determine  the 
sign  of  the  contribution.  Through  second  order  (in 
which  no  folded  diagrams  or  equivalent  complications 
arise)  the  sign  is  (-  1)"*,  where  h  is  the  number  of 
hole  lines  and  l  is  the  number  of  loops.  For  the  pur¬ 
pose  of  loop  counting,  double-  and  higher -excitation 
diamond  vertices  must  be  split  into  two  and  more  (re¬ 
spectively)  single  particle-single  hole  vertices,  with 
the  pairing  of  open  lines  corresponding  to  the  second- 
quantized  definitions;  for  example,  a  with  i  and  t>  with 
j  in  Eq.  (49).  A  factor  of  j  must  be  included  for  each 
pair  of  equivalent  lines.  Note  that  the  open  lines  of 
the  model  states  are  never  equivalent.  The  factors  of 
5  in  Eq.  (34)  are  in  addition  to  those  above. 

B.  Second-order  theory 

The  second-order  effective  Hamiltonian  is  given  as 

<V  =  5  <4>„  i  (K«  +  K  )FN<t>s )  .  (58) 

The  numerators,  which  comprise  most  of  the  informa¬ 
tion  explicitly  illustrated  in  a  diagram,  are  identical 
for  these  two  terms.  However,  in  second  order,  we  en¬ 
counter  the  first  summations  and  denominators  of  the 
resolvent.  We  may  rewrite  Eq.  (58) 


<V  =  \  L  (*JV'N4>i><4>,iFN>tV 

*  U*r> 

(59) 

where  E J  is  defined  by 

tf°4>r  =  £?4>t  , 

(60) 

and  <f>j  is  an  external  (i.e. ,  nonmodel  or  complementary) 
state  determinant.  We  find  that  if  4>a  is  a  single  exci¬ 
tation,  the  interactions  may  contribute  if  ■!>,  is  a  single, 
double,  or  triple  excitation.  If  <t>a  is  doubly  excited. 

may  even  be  a  quadruple  excitation.  Since  E ^  and 
El  both  contain  £?F,  the  denominators  are  functions  of 
only  those  orbital  energies  c  which  correspond  to  labeled 
particle  and  hole  lines  in  the  diagram.  For  example, 
the  numerator  of  the  matrix  element  between  two  sincly 


FIG.  1.  First-order  ettectivo  Hamil¬ 
tonian  matrix  element  diagrams. 

Lines  above  toelow)  the  horizontal 
midline  ot  the  diagrams  repres 
particles  ihotesl.  The  arret'  :  :■»  - 
ruled  for  elaritv,  are  redundant  -  .t .' 
this  notation. 
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excited  determinants  involving  a  singly  excited  inter¬ 
mediate  state  is 

E'wiv  ♦s><*siv«  (o>i\\ci)(cj\\bk) . 

C.» 

(61) 

The  prime  is  a  reminder  that  model  states  are  not  to  be 
included  in  the  sum.  The  "denominator"  factor  [which 
multiplies  Eq.  (61)1  is  given  by 

I  { [“- 

+  [<«*-  «/)  "(Z  “  Z  fe)j  |  -  (62) 


where  the  sums  over  h  and  p  designate  all  hole  and 
particle  lines  which  would  be  intersected  by  a  vertical 
line  drawn  between  the  two  interaction  vertices.  In 
the  case  of  a  singly  excited  intermediate  state,  these 
sums  consist  of  a  single  orbital  index  each. 

The  algebraic  expressions  for  the  second-order  ef¬ 
fective  Hamiltonian  matrix  elements  are  given  in  Table 
IV  and  their  diagrammatic  forms  are  shown  in  Fig.  2. 
To  save  space,  we  have  indicated  roughly  half  of  the 
required  diagrams  as  particle-hole  reversals  (phr)  of 
the  illustrated  diagrams.  From  top  to  bottom,  the  sub¬ 
divisions  correspond  to  the  HF  IHF,  HF ISE,  SE  ISE, 


TABLE  IV.  Second-order  effective  Hamiltonian  matrix  elements. 


<*0I  VNR|jvN  *„)  =  i  £  - 


{e/W  mnr 


*t<m  Dim 


i  <4\,l  V^o-t-Ra)  Vv*J)  =  •'  [  2  (cellmn)  {nm  lAe)  Km+ri -e-c)*'  +(m*n  —<?—*>*')  -  {phr)- 

(iXl 


1  I  VN  +Ra)VN  *J)  =  V  2  (am II re)  (me lie*)  Ua  +  m  -i  -e)‘‘  +  (c  +m  -*  -  el* 


+  •[  i  6,,  2  (amll/e)  (e/ll  me)  [(m  +  a -e  +  (m  +c -e  i  2  'ce  ®  nm')  <nm  *Dm<  I 

1  ttm 

-i  S  (*alle/)  (felltc)  ((*+a -e  -/)*’  +(f+c -e-/)’1)  +  [  2  (cellaro)  (feel!  mi)  ((c  +  m -e -c)"1 

0f 

ItfM  l«lml 

*(i  +  »>  — *  — eC'l!  •*  {phrs}  +  1  2  (*m  lice)  (ae  Him)  [D£,"'  +  i  Tj  — 

sm  *fmn  Urrm 


HFI  DE: 


J  (*0 1  VN  (Ho+Raiv'1  **')  *  V  2  ((dc)el1(f*)m)  (m(cd1  lle(*f>)  ((I  +  m  -d  -e)"1  +(c  +m  -*  -  of1! 

+  •!  i  2  '  (W  lief)  <e/lled)  ((*  +1  -e  +(c  +d  — c  -O'1)} +  (phr) 

•  «/ 
l«/»!l 

SE  I  DE: 

t’N  (R“  +  VN  *Jf)  =  ]  V  2  (aUk)l\ie>  (e(AD  IldcMia +f -/ -e>‘ 1  +  (c +d -* -e)*1)  J  +tphr> 

+  2  (Wlallep  (e/llcrf)  Ka +  <*/> -e -/T1 +  (c  +d-c-/)''l  -  V  2  M  Sfr'  \arii  cd 

4/  * 

(tf/Uf  *«***• 

^  |(fe — i  — i*)"1 +  — a  — r)",|+  J  2  v  » (JtHv  !1  (erf)/ 

*  ((<i  v-/ -  (t/c) -.^T1 +  (krf) -M) -t*)*1]  +  i  £  «7»i  Wc)<"  ,U’/)r  (crf)m 

1  (tfcWlml 

■  j(a  +  m -Wc* -c)*1  +  ((cd)  +m —  i*/))*'l  J  -{plirs-  +  |  -  ;  ^  it!  </'  >rf  .  (ed)i 

>  IDff1  ♦(f+(cd)-e-rr,)i '  £  t«  >>»  ■•»  *'■ 

v  I^TmL ♦(»'  +lcrf) -e  -O'1)  +  1  S,<4c)  2  '(cd)r  4Wli»'  im  ■  tf*lc  IDiJJJ’l  +ti -m -i/H -ef'! 


-  1  2  \cdll  mi/*)'  ((*f)m  Hail  +  (i  * m  -(*/>- n )*! I  |  -(phrs, 
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TABLE  IV  ( Continued ) 

DE  I  DE: 

t  |  *»r  *  |  J  23  il  (i/(e )  (\kl)e ilcdlKa  *b  -  (lit  -i’)'1  +  (c  *d-  Uri)  -c)"1 1 +iphr| 

wwni 

*i  6<jni<»>  S'  iiablm  II  (,/),')  (Well  (ed)m  l((ab>  *m _(,;) -e)'1  *Ucd)  *  m  -  W -ef'\ 

em 

+  |  J  6)tu  23  \ab\lef)  (e/Wcd)  lla+b -e -f)~' +(c  T-d-e -n~'\ 
iS/iil 


-J  5(ji)U»)  23*  (UiHdc)  II  (ab)ni )  ( (cd)m  il  MHbal  >  U»i  +  (ab)  -  (i/>  -  Wc»)*‘  +(m  *c(d)-(kl)~  Wall''  I 

(ttaXtfcXji )»"} 

+  i  6<>(,u*>  S'  a*  lleidci)  (e(kl)ll iijHcd)')Ha  +  h  —  c  —  irfc))*1 » (led)  +  (ij)  —  (hi)  —  r)"!l 


+  J  E  \(kL)(ab)\\e(i)))  (e(ba)\lcd)U(abl  +  (kl)  -Ui)-e)’'  <-lc  +d -e  ~(ba»-1]^  MPhrs? 

(<t«>rtil 

+  I  J  5'jw  o^hjci  E  (e/H  (a*)»i>  Ucdlm  tle/MUa*  •  *»<  —«  -/•*'  +  Ucd)  +  m -e  -/)** ! 

'  <?/ni 

-2  aj(»i  E  (cldc)  II  (ai)m  'l  ( (cd)infl  c(6a))  |((ai)  +  »;  —  e  —  (dc))"1  +  Kcd)  +  w  -  e  -  (ball’1 1 

em 

+  3  Stum,)  Knn,c>  E  ef 11  WtZH6a)>  (( cdHii )  lle/MUfta)  +  W -e  -/)*'  +  «cd)  +(;<)  -e  -/)''! 

ef 

+  £  5(i/)„n  23  U’WcMI  (aiMffeO  Ucd)(ji)ll  e(6a)>  IUa6)  +  (/*>-<•-  Wcl)*1  +  Kerf)  +  <»)  -e  -  (*a»)'1! 

e 

+  £  a<i»<»n  23  '.(cdHii)  H  Ilklm)  ((dc)mllab)  ((»«+«<>) -(yi)  -  <cd)l''  +  (ih  +  We)  -  A -a)‘l I 

m 

-4  23  <(»7>w  H  (ab)(lk))  Km  +  (ba)  -c  -rf)'1  *(m  *  (ij)  -Uk>  ~  (a*))*M 


-i  6wui>  E  (l«.i)?H  (m(dc)  II (aile)  ((m  +  tffr )  —  («/)  -W1  ■*•(»»)  +  (dc)  -e  -  lailfM  J  -\phrs} 

em 


s  *(„,(«,  *<,,<,><«>  E  II (6a>e-><e-<rfc» II  ^ 


efm*  De£ 

Mortal  E  (ISalellmaXa.alleldcIMU^^'H-D^^Vi  4»<  E  -fc 

f  emu  m 

£  a<«*i<<*>  E(0IHa6)e>(^WO«i*>  \  +{phrs} 


+{lf  4(|J)  's  a  model  SE,  Include 

-i((/im*a)llcd><(a4)«*)IU;)lUi»a)  +  («)  -c  -rf)'1  +(i  +  .i  -  Uk)  -  (a*))'1!}  *{phr[ 

+{H  *(*  is  a  model  DE  and  *J{  Is  not,  include  J (ij  life!  Kerf  II ab  'Kfe+f  -  i  - ))"'  +(e  +  rf  -  a  -i>f1!1l-,-{phr} 

+  If  *(S?'({J!  is  not  a  model  DE  and  ■*>( ?/>(*/>  is>  include  5  <<i6llcd'KM  IU.iMla  +  6  -c -rf)‘l  +  (i  »  I -h -D"1] 

*  1  < ij  II ah )  < cd II « ) (D?fl  +C,T'l 


HF  IDE,  SE  IDE,  and  DE  IDE  blocks,  TheSEIHF, 

DE  IHF,  and  DE  ISE  diagrams  could  be  obtained  by 
time  reversal  of  the  HFISE,  HFIDE,  andSEIDE  dia¬ 
grams,  respectively,  but  they  are  not  required  because 
of  the  Hermiticity  of  X,tt.  Except  for  the  denominator 
factors  (which  are  not  shown  in  the  figure)  and  the  ground 
state  term  on  the  diagonal,  the  SE  ISE  diagrams  are 
equivalent  to  those  of  Paldus  and  filzek.  Each  inter¬ 
mediate  state  (in  Fig.  2)  has  been  labeled  as  singly 
(S),  doubly  (D),  triply  (T),  or  quadruply  (Q)  excited, 
relative  to  <f>o. 


In  the  table,  {phrs}  represents  the  phr  contributions 
analogous  to  the  terms  in  the  immediately  preceding 
brackets.  Except  that  the  phr  terms  are  grouped  to¬ 
gether,  the  order  in  the  table  corresponds  to  that  in 
Fig.  2.  The  algebraic  expressions  for  the  phr  terms 
are  obtained  by  exchanging  (relabel  mg)  a  and  i.  b  and 
j,  etc.  in  the  illustrated  terms.  Sign  changes,  when 
required  (due  diagrammatically  to  a  change  in  the  parity- 
factor  with  particle  hole  reversal  when  the  numbers 
of  particles  and  holes  are  not  both  even  or  both  odd), 
have  been  explicitly  indicated.  Since  the  type  of  inter- 
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FIG.  2.  Second-order  effective  Hamiltonian  matrix  element 
diagrams. 


mediate  state  is  not  always  obvious  in  the  algebraic 
form,  the  superscripts  and  subscripts  of  the  terms  to 
be  omitted  from  the  primed  sums  when  they  correspond 
to  model  states  have  been  enclosed  in  curly  brackets 
below  the  summation  sign.  The  level  of  excitation  is 
determined  by  the  number  of  labels  provided. 


Note  that  the  parenthetic  notations  sometimes  appear 
without  delta  functions.  In  these  cases,  they  indicate 
multiple  terms  which  must  all  be  included.  Terms  con¬ 
taining  parenthetic  notations  are  indicated  in  the  figures 
by  singie  diagrams,  since  with  unlabeled  lines  their 
skeletons  are  identical.  The  negative  double  excitation 
energy  notation 

=  +e,  (63) 

has  been  employed  within  the  table,  and  orbital  energies 
have  been  indicated  by  their  subscripts  alone. 

We  see  that  the  second-order  (closed)  ground  state 
correlation  contribution  appears  along  the  diagonal. 
These  terms  are  considered  linked,  since  noninteracting 
open  lines  merely  represent  delta  functions  which,  in 
this  case,  indicate  that  the  diagram  contributes  only  to 
diagonal  elements.  They  arise  from  triply  excited 
intermediates  for  singly  excited  determinants,  and  from 
quadruply  excited  states  for  doubly  excited  model  states. 
Since  each  intermediate  (including  the  open  line  parts) 
is  doubly  excited  relative  to  the  model  state,  the  de¬ 
nominator  dependence  on  the  model  state  is  removed  and 
the  denominators  are  identical  to  those  of  the  ground 
state  terms.  Assuming  that  the  model  space  does  not 
contain  any  triply  or  quadruply  excited  states,  the 
value  of  this  term  is  exactly  equal  to  that  obtained  from 
ordinary  MBPT,  except  that,  in  the  HF  !HF  matrix  ele¬ 
ment,  the  contributions  from  any  doubly  excited  deter¬ 
minants  included  in  the  model  space  must  be  omitted. 
These  omitted  terms  correspond  to  the  couplings  (i.e. , 
off-diagonal  effective  Hamiltonian  matrix  element)  be¬ 
tween  $0  and  the  doubly  excited  model  states,  which  are 
explicitly  included.  The  identical  terms  contributing  to 
the  excited  state  diagonal  elements  are  not  omitted. 
Subtraction  of  the  full  ground  state  diagrams  from  the 
diagonal  gives  the  effective  excitation  energy  matrix 
relative  to  the  one-dimensional  MBPT  ground  state 
energy. 

In  the  second  order,  the  diagrammatic  differences 
between  general  and  complete  model  space  implementa¬ 
tions  appear.  We  find  that  because  the  general  methods 
allow  the  Q  space  to  include  states  that  are  related  to 
(i.e. ,  consist  of  particle  and  hole  labels  in  common 
with)  the  selected  P  space  states,  certain  disconnected 
“irreducible”30  diagrams  must  be  included.  These  dia¬ 
grams  consist  entirely  of  open  lines,  and  thus  they  are 
trivial  to  calculate  because  no  summations  are  required. 
In  addition,  many  of  them  cancel  each  other.  These  dis¬ 
connected  diagrams  are  those  in  the  last  two  rows  of 
Fig.  2  for  DE  IDE  matrix  elements.  The  last  diagram, 
with  <t0  as  the  intermediate  state  labeled  X,  is  not 
used  when  the  model  space  includes  $0.  However,  if 
4>o  is  included  in  Q  rather  than  P,  this  diagram  is  con¬ 
sidered  and  is  found  to  cancel  the  immediately  pre¬ 
ceding  quadruple  intermediate  diagram  exactly,  so  that 
neither  diagram  need  be  calculated.  In  fact,  close 
examination  reveals  that  none  of  these  disconnected  dia¬ 
grams  would  remain  if  the  complication  of  excluding  the 
model  states  as  intermediates  did  not  arise.  We  find 
for  the  case  that  the  model  space  has  no  singly  excited 
states  that  the  single  and  triple  intermediate  diagrams 
of  the  penultimate  row  cancel  each  other.  However,  if 
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the  model  space  does  include  singly  excited  deter¬ 
minants,  the  corresponding  terms  must  be  omitted  from 
the  S  diagram.  Assuming  that  the  model  space  does  not 
include  triply  excited  states,  the  related  terms  from  the 
T  diagram  do  appear  and  must  therefore  be  calculated, 
since  they  are  not  cancelled.  Thus,  in  Table  IV  the 
sixth  term  from  the  end  (counting  phrs,  i.  e. ,  the  first 
“conditional"  term)  gives  the  net  contribution  of  S  and 
T  in  the  presence  of  model  singles. 

Similarly,  the  last  two  D  diagrams  need  only  be  con¬ 
sidered  for  special  cases  involving  model  intermediates. 
The  first  is  usually  cancelled  by  its  phr  diagram.  How¬ 
ever,  there  may  be  instances  in  which  the  cancellation 
is  not  permitted.  For  example,  (assuming  that  and 
are  model  states)  if  is  a  model  determinant  but 
is  not,  the  term  with  as  intermediate  must  be 
calculated  while  the  term  with  is  omitted.  The  re¬ 
maining  D  form  represents  sixteen  individual  terms, 
which  cancel  with  each  other  unless,  for  example, 
is  a  model  state  and  is  not. 

The  last  six  diagrams  are  not  required  if  a  "com¬ 
plete”  model  space  is  specified,  i.  e. ,  if  P  contains 
all  possible  excitations  among  a  group  of  “active"  or¬ 
bitals  { i,  j  •  •  •  -  a,  b  •  •  • }.  In  this  event,  the  dia¬ 
grammatic  expansion  has  been  described  as  “fully 
linked.  "  Unfortunately,  an  undue  stigma  has  been  as¬ 
sociated  with  the  disconnected  open  diagrams  which 
must  be  included  for  the  general  ("incomplete”)  model 
space.  These  terms  are  necessary  and  are  size  exten~ 
sive.  Hose  and  Kaldor  have  discussed  the  size  ex- 
tensivity  of  similar  terms  which  arise  in  their  meth¬ 
od.30-38 

Although  the  present  diagrams  appear  to  be  quite  dif¬ 
ferent  because  we  have  chosen  to  show  the  explicit  dia¬ 
grams  for  the  chosen  model  space  and  to  include  all  of 
the  open  lines  in  them,  the  diagrams  we  require  are 
analogous  to  the  pure  two-electron  interaction  forms 
given  by  HK,  30,38  except  for  our  denominator  combina¬ 
tions.  The  complete  equivalence  through  third  order, 
achieved  by  constructing  Hermitean  combinations  of 
their  terms,  may  be  easily  verified.  If  one  does  not 
make  Hermitean  combinations,  the  non-Hermitean  for¬ 
mulation  requires  that  the  entire  effective  Hamiltonian 
matrix  be  calculated,  rather  than  its  upper  (or  lower) 
triangular  part.  Additionally,  we  find  that  the  cancella¬ 
tions  described  above  for  our  disconnected  diagrams 
do  not  occur  for  the  non-Hermitean  form. 

It  appears  that  we  have  more  diagrams  than  HK,  but 
this  is  because  we  explicitly  show  the  various  labelings 
of  the  open  diagrams  that  must  all  be  included  in  either 
method.  Also,  since  our  particle  and  hole  states  are 
not  redefined,  unlike  the  HK  t'ansition  method,  it  is 
easier  to  recognize  the  model  states,  which  in  either 
method  must  be  omitted  as  intermediates.  For  ex¬ 
ample,  the  similar  appearance  of  and  of  a  quadruple 
excitation  from  <t>v  when  represented  relative  to  a 
double  excitation  from  can  be  confusing.  Further¬ 
more,  with  our  convention,  a  given  model  state  has  one 
definition  rather  than  a  set  of  definitions  of  the  dimen¬ 
sion  of  the  model  space.  This  should  simplify  program¬ 
ming. 


The  provision  for  deleting  the  possibly  large  numbers 
of  more  highly  excited  model  states  from  P  is  an  attrac¬ 
tive  practical  feature  of  general  model  space  methods. 
This  greatly  reduces  the  number  of  effective  Hamiltonian 
matrix  elements  to  be  evaluated,  and  also  the  number 
of  potential  intruder  states.  The  method  is  not  restricted 
to  the  single  and  double  excitation  model  space  illus¬ 
trated  in  this  paper,  and  could  be  applied  to  other  prob¬ 
lems. 
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Equations  for  the  first-  and  second-derivative  coupling  matrix  elements  between  the 
adiabatic  electronic  wave  functions  of  a  canonical  Van  Vleck  approach  are  derived  for  the 
case  of  a  single  nuclear  coordinate.  Expressions  for  their  order-by-order  evaluation, 
analogous  to  the  multidiT-  nsional  perturbative  expansion  of  the  effective  Hamiltonian 
(given  elsewhere)  are  presented.  Diagrammatic  notations  useful  for  enumerating  certain 
components  of  the  derivatives  are  introduced,  and  the  diagrams  required  through  first 
order  are  shown. 


I.  INTRODUCTION  j 

The  Bom-Oppenheimer  adiabatic  approxima*  I 

tion1  is  one  of  the  most  fundamental  and  widely  I 
used  approximations  in  the  theory  of  molecular 
structure.  When  the  electronic  wave  function  is  ( 
slowly  varying  with  respect  to  nuclear  displace-  j 

ments,  solution  of  the  electronic  Schrodinger  equa-  : 
tion  leads  to  accurate  calculations  of  the  properties  j 
of  molecules.  Adiabatic  potential  energy  curves 
and  surfaces  are  used  in  the  interpretation  of  spec¬ 
troscopic  data,  the  prediction  of  thermochemical 
quantities  (heats  of  reaction  and  formation  and  ac¬ 
tivation  energies),  and  in  the  approximate  descrip¬ 
tion  of  the  dynamics  of  molecular  collisions  (at  en¬ 
ergies  of  chemical  interest)  usually  involving  a  sin¬ 
gle  electronic  state.  Nonadiabatic  treatments 
(direct  solution  of  the  unseparated  Schrodinger 
equation)  have  also  been  pursued.2  The  Bom-  ; 
Oppenheimer  adiabatic  representation  can  still  be 
used  to  accurately  describe  dynamical  phenomena, 
such  as  radiationless  transitions  between  electronic 
states  of  bound  molecules  or  electronic  transitions 
in  molecular  collisions,  by  consistently  including 
the  nonadiabatic  coupling  terms  in  the  process  of 
solving  the  nuclear  Schrodinger  equation.3  The 
most  important  coupling  terms  for  many  systems 
(and  those  generally  most  difficult  to  obtain)  are 
the  matrix  elements  of  the  nuclear  derivative  opera¬ 
tor.4 

Recently  there  has  been  increasing  interest  in  the 
calculation  of  nondiabatic  couplings  for  ab  initio 
electronic  states.5-11  Analytical  methods  of 
evaluation  must,  of  course,  be  derived  specifically 
for  each  electronic-structure  method.  Those  ad¬ 
dressed  thus  far  are  the  Hartree-Fock  self-consis¬ 
tent-field  (SCF)  and  configuration  interaction  (Cl)  , 
methods.  Investigations  of  the  utility  of  the 


Hellmann-Feynman  theorem5  or  the  Sidis  formula6 
(both  valid  for  exact  states)  in  the  case  of  practical 
(inexact)  wave  functions  have  been  made.  Numeri¬ 
cal  methods  based  on  finite-difference  formulas,7-8 
although  applicable  to  any  type  of  wave  functions, 
need  repetitive  solutions  of  the  electronic  problem, 
a  nontrivial  requirement  when  highly  correlating 
methods  are  employed.  The  approach  of  fitting  a 
polynomial  followed  by  analytical  evaluation  of  the 
derivatives8-9  becomes  cumbersome  for  large  Cl 
wave  functions  and/or  for  systems  with  high  coor¬ 
dinate  dimensionality. 

The  present  paper  derives  analytical  expressions 
for  the  nuclear  derivatives  between  electronic  wave 
functions  employed  in  multidimensional  diagram¬ 
matic  perturbation  theory.  As  in  that  theory,  ex¬ 
plicit  calculation  of  the  wave  functions  is  unneces¬ 
sary.  The  approach  uses  a  formulation  for  the 
first  derivatives  that  avoids  numerical  differentia¬ 
tion  procedures  and  provides  an  order-by-order 
prescription  for  obtaining  increasingly  accurate 
j  values.  The  theory  incorporates  the  effects  of  con¬ 
figurations  omitted  from  the  model  space  and  is,  in 
principle,  exact. 

For  comparison,  the  expressions  for  the  matrix 
elements  of  the  first-  and  second-derivative  opera- 
1  tors  using  Cl  wave  functions  are  given  in  Sec.  II. 

In  Secs.  Ill  and  IV  the  nonadiabatic  coupling  ma- 
j  trix  elements  corresponding  to  the  perturbative  ex- 
I  pansion  of  the  electronic  wave  function  are  de- 
!  rived.  In  Sect  V  the  evaluation  of  these  matrix  ele- 
j  ments  at  the  orbital  level  for  a  specific  illustrative 
model  space  is  presented  in  the  conventional  and 
1  second-quantized  notations,  and  diagrammatic  no¬ 
tations  for  the  required  elementary'  components  are 
introduced.  The  derivation  of  an  exact  Hellmann- 
1  Feynman-like  relation  is  included  as  an  Appendix. 


II.  DERIVATIVE  MATRIX  ELEMENTS 
IN  THE  Cl  FORMALISM 

Solution  of  the  full  Schrodinger  equation  can  be 
accomplished  by  separating  the  Hamiltonian  into 
nuclear  and  electronic  pans 

H=ffnu+tfd  (1) 

and  assuming  a  sum-of-products  form  of  the  total 
wave  function 

*;,0,=  £  <(/;;£)*£(&).  !  (2) 

r 

We  use  the  subscripts  j  a,f3,y. . .  )  to  indicate  the 
members  of  the  set  of  states  of  interest.  The 
Bom-Oppenheimer  functions  ( 'Pj! )  are  solutions 
of  the  electronic  Schrodinger  equation 

(f/d~£d)'P'1=:0.  (3) 

In  principle,  a  complete  set  of  ^’s  is  required,  but 
for  practicality  it  is  assumed  that  a  reasonably 
small  number  of  electronic  states  { t by;y=a,0. . .  J 
will  form  a  sufficient  expansion.  In  the  Cl 
method,  the  electronic  wave  functions  may  be  ex¬ 
pressed  in  a  basis  of  determinants  (or  of  configura¬ 
tion  state  functions) 

,  (4) 

<r 

where  j  cx,r. . .  ]  are  labels  for  a  complete  expan¬ 
sion.  The  coefficients  (  Cay  |  are  solutions  to  the 
secular  problem 

HC  =  CE'[ ,  (5)  ; 

where 

H  =  {  (6) 

and  |  <t>„  j  are  orthonormalized  functions,  such  as  j 
Slater  determinants  formed  from  Hartree-Fock  or-  J 
bitals.  An  explicit  dependence  on  the  nuclear  coor¬ 
dinates  R.  arises  because  these  coefficients  as  well 
as  the  expansion  coefficients  (not  shown)  of  the  d>’s 
in  terms  of  the  primitive  atomic  orbital  (AO)  basis 
are  different  for  each  nuclear  arrangement.  An 
implicit  g.  dependence  of  4  results  from  the 
nuclear-centered  nature  of  the  AO  basis.  The  best 
solution  for  a  given  AO  basis  is  obtained  when  the 
“full  Cl  problem”  is  solved;  that  is,  the  sum  in  Eq. 
(4)  (and,  thus,  the  dimension  of  H)  is  over  all  pos-  ! 
sible  determinants  (excitations)  within  the  basis.  i 
However,  in  practice,  the  Cl  problem  is  almost  al¬ 
ways  truncated. 

For  simplicity,  we  examine  the  form  of  the  ma¬ 
trix  elements  of  the  nuclear  Hamiltonian  for  the  ! 


case  of  a  diatomic  molecule.  These  couplings  in¬ 
volve  derivatives  of  the  electronic  wave  functions 
with  respect  to  the  single  intemuclear  coordinate 
R,  (^l  B'Pg/BR  >,  and  ('Pa  j  B2'V$R2),  where  we 
have  dropped  the  superscript  on  )Pri.  In  the  Cl 
framework  the  first  derivative  is  given  by 


('I'° 


a 

a  r 


'<>0 


+2cl(<t>a 


(7) 

where  we  use  the  convention  that  the  derivative 
operator  may  not  operate  beyond  the  ket  of  an  in¬ 
tegral. 

It  is  well  known  that  the  coefficient  derivative  in 
the  first  term  become  large  in  regions  of  avoided 
crossings  (between  states  of  the  same  symmetry). 
This  might  cause  difficulties  for  a  numerical 
evaluation  scheme.  If  so,  it  is  expected  that  the 
matrix  elements  of  the  Hamiltonian  would  be  more 
well  behaved  than  the  coefficients,  so  that  numeri¬ 
cal  methods  based  on  an  exact  Hellmann- 
Feynman-like  relation  (see  the  Appendix) 


9  OyT 


JL  R 

BR  aT 


(8) 


could  be  used.  In  this  notation,  the  derivative 
operator  must  stop  at  the  right  parenthesis,  and  Eg 
is  a  specific  element  of  £d.  Although  Eq.  (8)  is 
valid  for  any  wave  function  satisfying  Eq.  (5),  this 
approach  is  cumbersome  since  obtaining  accurate 
energies  requires  the  dimension  of  ft.  to  be  large; 
e.g.,  the  basis  of  Eq.  (4)  might  consist  of  a  given 
set  of  determinants  and  all  single  and  double  exci¬ 
tations  from  that  set. 

The  integral  in  the  second  term  is  zero  for  a^r 
if  <t>„  and  <t>r  differ  by  more  than  one  orbital,  and 
for  cr  =  r  if  the  <t>’s  are  real.  Some  treatments  have 
assumed  that  these  terms  vary  slowly  with  R  and 
have  neglected  them.  In  fact,  however,  they'  can  be 
of  a  magnitude  comparable  to  the  first  term.9 
When  the  4>*s  are  chosen  to  be  Hartree-Fock  deter¬ 
minants,  the  derivatives  in  the  second  term  can  be 
expressed  in  terms  of  quantities  used  in  and  ob¬ 
tained  from  the  coupled  perturbed  Hartree-Fock 
(CPHF)  method,11,1'  which  is  commonly  applied 
for  the  determination  of  energy  gradients.  Alter¬ 
natively,  they  have  been  obtained  by  numerical9  or. 


4 


in  some  cases,  analytical711  means.  For  the  pur¬ 
poses  of  the  present  paper,  we  shall  assume  that 
the  derivatives  of  the  form  in  the  second  term  are 
available  and  refer  to  them  as  standard  or  configu¬ 


ration  derivatives. 

The  second  derivatives  are  more  complex. 

Direct  differentiation  of  Eq.  (4)  leads  to  the  expres¬ 
sion 


An  alternate  formulation  is 


(10) 

where  the  first  term  would  be  obtained  numerical¬ 
ly.  If  the  full  Cl  problem  is  solved,  the  second 
term  can  be  expanded  by  insertion  of  the  resolution 
of  the  identity 

i— 2  i^x^l 

o 

=  2  l  V<^l  +  2  IV<*Vl  (11) 

r  <* 

(where  the  subscripts  { fi,v. . .  }  indicate  members 
of  the  complement  to  (  a,0 ...  1 )  to  obtain  an  ex¬ 
pression  entirely  in  terms  of  first-derivative  matrix 
elements14: 


(12) 


Otherwise  this  formula  is  only  approximate.  Note 
that  it  requires  additional  first  derivatives 

((•’•Im*-)]- 

i.e.,  couplings  between  the  basis  states  of  interest 
and  those  “not  of  interest.”  The  effect  of  trunca¬ 
tion  of  this  sum  might  be  particularly  severe  be¬ 
cause  the  first  term  of  Eq.  (10)  includes  the  second. 
Thus,  for  an  approximation  to  Eo.  (12),  the  intend¬ 
ed  cancellation  might  be  incomplete.  This  would 
cause  spurious  contributions  from  the  first  term  to 
remain. 


III.  DERIVATIVE  MATRIX  ELEMENTS  IN  THE 

CANONICAL  VAN  VLECK  FORMALISM 

In  the  past  few  years  methods  for  the  determina-  ' 
tion  of  potential  energy  surfaces  have  envolved, 
such  as  direct  Cl,15  unitary  group  Cl,16  single-state 
many-body  perturbation  theory  (MBPT),n,ls  mul- 
I  ^configuration  SCF  (MCSCF),19  and 
'  MCSCF  +  Cl.19  For  many  applications  MBPT 
has  several  advantages  over  the  Cl  approaches.  j 
However,  a  perturbative  treatment  based  on  a  sin¬ 
gle  (SCF)  determinant  is  restricted  to  the  lowest  \ 
state  of  a  given  symmetry.  It  may  even  be  expect-  j  j 
ed  to  encounter  difficulties  for  that  state  in  regions 
of  avoided  crossings,  where  more  than  one  configu¬ 
ration  is  important.  In  addition,  this  method  is 
not  well  suited  for  many  applications  involving 
open  shells  or  multiple  bonds.  One  solution  to 
these  problems  is  the  canonical  Van  Vleck  ap¬ 
proach  to  quasidegenerate  perturbation  theory.10 

When  the  canonical  Van  Vleck  formalism  is  im¬ 
plemented  via  a  diagrammatic  expansion,2'  the  en¬ 
ergies  obtained  from  approximate  solutions  of  Eq. 

(3)  will  be  size  extensive,22-23  unlike  energies  of 
non-full-CI  calculations.  In  addition,  perturbation 
theory  provides  a  means  of  decreasing  the  size  of 
the  Hamiltonian  matrix,  which  must  be  diagonal- 
)  ized,  and  a  prescription  for  obtaining  increasingly 
.  accurate  (and,  in  principle,  exact)  results. 

In  this  formalism,  the  (exact)  correlated  electron- 
■  ic  state  of  interest  are  given  by 

/ 

Vr=2tf(£)<i>8 (G&tfVfi) .  (13) 

s 

1 

where  the  basis  for  expansion  has  been  restricted  to 
i  a  selected  set  of  states  j  4>j  }  corresponding  to  the 
states  of  interest  and  referred  to  as  the  model 
space.  The  £  dependence  will  not  be  indicated  in 
subsequent  equations.  The  j  D»r  )  are  the  eigen¬ 
vectors  of  the  effective  Hamiltonian: 

^rr£=££el .  (14) 

and  (  U  |  <t>8)  |  are  the  perturbed  (improved)  model 
functions.  U  is  the  unitary  decoupling  operator 


I 


25 


which  transforms  the  Hamiltonian  operator  H  to 
block  diagonal  form  £?.  In  principle,  the  exact 
eigenvalues  will  be  obtained  because  the  effects  of 
the  complementary  states  will  be  contained  in  U. 

In  practice,  U  will  be  determined  such  that  the 
transformed  Hamiltonian  operator 

jr=UfH''U  (15) 

is  block  diagonal  through  a  certain  order.  Then 
the  effective  Hamiltonian  JK.rr.  formed  of  matrix 
elements  of  JK,  |  <<t>„  |  ),  will  be  diagnonal- 

ized. 

Actually,  an  order-by-order  prescription  for  Jf*.rr 
is  given,  and  U  is  not  explicity  calculated.  Howev¬ 
er,  for  comparative  and  derivational  purposes,  the 
formulation  of  the  wave  functions  may  be  exam¬ 
ined.  Inserting  the  resolution  of  the  identity  in  the 
form 


1=2  I 

a 

=  2  lV<*r!  +  2 


into  Eq.  (13),  we  can  see  that  Vr  may  be  written 

a.  6 

Comparing  Eqs.  (4)  and  (17)  we  find  the  relation 
between  the  full  Cl  coefficients  and  the  exact 
(infinite-order)  effective  Hamiltonian  eigenvectors: 

.  <11 

6 

Starting  from  Eq.  (13),  we  find  (since  U  is  uni¬ 
tary)  that  the  first  derivatives  are  given  by 


J  2  ^0^  ~  2 Day  ^r0  +  22^ar^r6 


('*'a }  as  **)“?1>-ra R  Dr* 


+  2D«y?y*D»  • 
r.S 


re 


& '  can  be  thought  of  as  matrix  elements  in  the  ori¬ 
ginal  d>  basis  of  a  transformed  derivative  operator, 
or  as  derivative  matrix  elements  in  the  basis  of  per¬ 
turbed  model  functions  l/f J>- 
With  the  reduced  dimensions  of  this  approach  it 
should  be  practical  to  calculate  the  first  term  of 
Eq.  (19)  by  the  analog  to  Eq.  (8): 

2Z>«r3*i>*s*(£*~£«)  Dht 


Again,  this  expression  is  an  exact  equality  (at  any 
order)  because  of  Eq.  (14).  However,  the  exact 
couplings  will  only  be  obtained  [from  Eq.  (19),  op¬ 
tionally  incorporating  Eq.  (21)]  when  the  infinite- 
order  effective  Hamiltonian  is  used.  Although  the 
summations  in  these  equations  run  over  only  the 
model-space  dimensions,  the  effort  of  evaluating 
the  couplings  is  hidden  in  the  determination  of 
The  expansion  of  this  term  will  be  examined  in  de¬ 
tail  in  Sec.  IV.  We  note  that  only  the  lower  tri¬ 
angular  parts  of  gj  and  dJ^.a/dR  need  be  calcu¬ 
lated  since  is  anti-Hermitian  and  is  Her- 
mitian.  I 

The  matrix  elements  of  the  second-derivative 
operator  obtained  by  direct  differentiation  are  I 


+2^^rfr£>W  * 


•^*(^^—2^).  '  (23) 

Equation  (22)  can  be  expressed  alternatively  [by  substituting  Eq.  (13)  into  Eq.  (10)]  as 

l'1’8! a*T'|/p)  =  a*  jRDar  TRDr0+%  lKDar 

-2^(2^^2^^W-  (24) 

r.6  t  ? 

Equation  (12)  could  be  used,  truncated  to  a  sum  over  y  only  (since  the  remaining  states  are  not  determined 


i,.1 


in  the  effective  Hamiltonian  method),  but  the  previously  stated  criticism  of  this  procedure  would  apply. 

The  second  derivative  of  the  coefficients  in  Eq.  (22)  can  be  alternatively  obtained  by  similar  techniques.  We 
write  it  as 


2D 


♦  jld  _j_ 

aYdR2re~M 


^DardRDr0 


(25) 


where  the  first  term  is  to  be  obtained  numerically.  Inserting  the  identity  into  the  second  term  and  substitut¬ 
ing  Eq.  (21)  gives 


a  r 


(■ Ee-EaV  '2Dlr 
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a 

i 

a  r** 
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D&g 
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+  2  < E6-EarlDay 

r.8,«.  f,ij 


D^Ep—Ef, 


'Dl 


(26) 


The  latter  process  can  also  be  applied  to  the  second  > 
term  of  Eq.  (24).  The  choice  between  Eqs.  (22) 
and  (24)  is  then  further  influenced  by  observing 
that  Eq.  (24)  requires  the  additional  matrix  ele¬ 
ments  |  ),  while  Eq.  (22)  calls  for  matrix 

elements.  These  off-diagonal  j£'  matrix  elements 
can  be  obtained  analogously  to  the  model-space  di¬ 
agonal  block  of  however,  there  will,  in  general, 
be  many  more  of  them.  The  perturbative  determi¬ 
nation  of  £  and  g."  will  be  given  in  Sec.  IV. 


an  is  given  by20 

c*A  0)  c  r»0 
yb  - y  y 

^={<Pr\V<t>6)  , 

^4,=  j<<I>y|[K,Gn,]<l>6>  ,  (29) 

where 

G(t>,=0 , 


IV.  PERTURBATIVE  EXPANSIONS 

The  first-derivative  matrix  elements  in  the  effec¬ 
tive  Hamiltonian  formalism  are  prescribed  by  Eq. 
(19).  The  coefficient  derivatives  (first  term)  can  be 
calculated  numerically  or  from  derivatives  of 
via  Eq.  (21).  In  the  latter  case,  an  ordered  evalua¬ 
tion  scheme  emerges.  In  the  canonical  Van  Vleck 
method  it  is  convenient  to  define  the  transforma- 


tion  U  by 

U=ea, 

(27) 

with  1 

G  =  -Gf . 

(28) 

Then,  through  third  order,  the  effective  Hamiltoni-  j 


G‘”=2f*?Kpr-/>rra?!  ,  (30) 

r 

G,2,  =  2i  H?[V,Gu,]Pr-/>r[K,G"’]R“  }  , 
r 


and 


Ht[=H°+V, 
ff°<Pr=E}<t>r , 
Pr=|<J»rX«J>J,|  , 

„0  -c 

Kr~Z  IE°-E° 

it  ^>t 


(31) 


It  should  be  noted  that  G  has  no  block  diagonal 
(D)  parts: 


Cd  =  X  l<X>yX<t>ri  C7<^5X<I>6|  +2l<V<<,VlC<t>*><‘fv:  *0-  (32) 

r.s  it.v 

Thus  G  is  completely  specified  by  its  off-diagonal  part  Gy: 

=  \+r>l+r\°%H+n\  +  lV«Vl<7*r><*rl  I  • 

r.M 


(33) 


It  is  then  straightforward  to  see  that  the  derivatives  of  JEOr  involve  only  derivatives  of  the  two-electron  in¬ 
tegrals  and  the  (zeroth-order)  orbital  energies: 

dR  76  ~^bR  r’ 

r*‘)K 

8  1  3  |~  [<»r!^X»<l|y0t,><»>|y»t> 

a/J  **  ~  2  a/?  ~  (£°-£°)(£°-£°)  +  (£g-£“)(£?-£“) 

_  -C.  (4>rl  ^ ><<»,!  K<DmX«I>4|  V<X>h) 

(£S-£2)(£2-£°)  +  (£j-£°)(£°-£°) 

These  derivatives  may  be  determined  by  the  CPHF  method  as  implemented  by  Pople  and  co-workers. 12 

£  may  also  be  examined  in  an  ordered  fashion.  When  V  is  expanded,  j£  may  be  given  by  a  Hausdorf 
expansion: 

*--(^||^+[^<;]+i([^<I]’0|+-h)'  ,33‘ 

where  the  derivative  operator  is  allowed  to  work  all  the  way  through  to  the  ket.  Then  we  find 


When  G  is  expanded  by  order  we  find  the  low  orders  of  to  be 


Thus,  by  multiplying  Eqs.  (34)  and  (37)  by  the  appropriate  eigenvector  components,  prescribed  by  Eqs. 
(19)— (21),  the  coupling  matrix  elements  can  be  obtained  to  an  order  consistent  with  that  of  the  effective 


'I 


Hamiltonian.  Explicit  contruction  of  the  wave  functions  is  unncessary. 

The  transition  from  the  notation  above,  involving  operator  derivatives,  to  the  basic  expressions  is  straight¬ 
forward,  although  fairly  complex.  However,  a  shortcut  is  available  for  those  terms  involving  the  diagonal 
block  components,  such  as  the  leading  terms  above.  Since  GD  =  0, 

-^-<<t>r|C<i>6>=0, 

~<4>„|G*„>«0,  ,38) 

and,  therefore. 


IIH*H 

((^-li0*^^-!0^*-)) 

These  equations  are  valid  for  each  order.  For  example,  at  First  order 
(  I  <*rl 

^  (El-El)  (£°-£°)  ‘  m 

The  integrals  {  <  <t>y  |  F^}  )  are  the  same  ones  required  for  &*$},  and  the  derivative  matrix  elements  are  a 
subset  of  the  standard  ones  used  in  Eq.  (7).  In  this  form,  considering  the  derivative  operator  as  a  one- 
electron  operator,  g/* 11  can  be  easily  diagrammed,  as  will  be  demonstrated  in  Sec.  V.  The  nested  derivatives 
involving  ( dG/dR  )D  may  be  expressed  similarly. 

However,  terms  with  (dG/dR  )*  must  be  handled  directly,  substituting  and  differentiating  the  components 
of  Eqs.  (30)  and  (31): 


-  -  ( 1^^ ) 

U*r  > 

,  \dR  /  ,3 


-t2  2- 


(£°-£°)  T 


(£j-£®’ 

<*m! 


,  JL  . 

<<M 

v%) 

<4»^ !  Fd>5> 

+  a  r 

(£°- 

-£°)  jj 

(El-El) 

,  M  • 

k' 

1  V 

) 

-4-X 

'  F<J>5> 

o  a. 

t*t 

1 

+  3J? 

;  (£?-£*) 

<4>r|  V%) 


*m>(^  ~4>^<<*>v|F<b4> 

(£°— £°  )(£g— 


Although  this  result  is  fairly  complicated,  the  terms  involving  the  two-electron  integral  and  energy  deriva¬ 
tives  are  the  same  ones  required  for  d%*M/dR  [in  Eq.  (34)].  The  remaining  terms  of  1:1  can  be  di¬ 
agrammed  similarly  to  the  first-order  terms,  although  there  are  quite  a  number  of  them. 

The  .£"  matrix  elements  called  for  in  Eq.  (22)  can  also  be  examined  in  an  ordered  fashion  by  converting 
from  U  to  G: 


IF’0  +!  sF’°  ■«  + •■•)**)• 


/ 


Then,  expanding  G  by  orders  and  collecting  terms  we  obtain 


3  R1 

A +2  (*' 

3  R 

JLg<2> 

3  Rz 

3  R 

Examples  of  the  evaluation  of  £'  and  in  terms  of  orbitals  are  given  in  Sec.  V. 


V.  diagrammatic  representation 

We  present  an  orbital  formulation  of  the  nonadi- 
abatic  couplings  between  electronic  states  obtained 
from  a  diagrammatic  implementation  of  canonical 
Van  Vleck  perturbation  theory,  illustrated  for  the 
model  space  chosen  in  Ref.  (21),  which  consists  of 
a  Hartree-Fock  determinant  plus  selected  single 
and  double  excitations.  The  effective  Hamiltonian 
matrix  elements  have  been  given  both  formally20 
and  diagrammatically.21  These  diagrams  are  useful 
for  showing  which  subsets  of  integrals  are  required 
for  evaluation  of  the  matrix  elements  and  (as  in  or¬ 
dinary  MBPT  calculations)  for  providing  size- 
extensive  approximations  and  efficient  algorithms 
for  their  evaluation.  The  derivatives  of  these  ma¬ 
trix  elements  required  for  the  evaluation  of  the 
coefficient  derivative  part  of  the  nonadiabatic  cou¬ 
plings  using  Eq.  (21)  will  be  obtained  via  CPHF,  as 
mentioned  earlier.  The  terms  in  £  can  be  di¬ 
agrammed  using  the  one-electron  nature  of  the 
derivative  operator  working  on  a  Slater  deter¬ 
minant.  Examination  of  the  derivatives  of  £ <0’, 
i.e.,  the  standard  configuration  derivatives 


(*r|^ 


in  the  conventional,  second-quantized,  and  di¬ 
agrammatic  forms  serves  as  a  convenient  introduc¬ 
tion  to  the  many-body  techniques  employed  herein. 

It  is  well  known  that  the  configuration  deriva¬ 
tives  in  Eq.  (36)  [or  Eq.  (7)]  are  zero  if  y  and  6  (or 
a  and  r)  are  identical  and  $  is  real,  or  if  <t>r  and 
<I>j  differ  by  more  than  one  orbital.  For  example, 
let  <t>r  and  06  be  singly  excited  determinants, 
which  may  be  expressed  in  terms  of  annihilation 
and  creation  operators  (indicated  by  orbital  labels, 
appropriately  daggered)  working  on  the  Hartree- 
Fock  determinant  <J>0: 

(4>y)=ati  ,  (44) 

\*s)~bfj\%) , 

where  a,  i,  b,  and  /,  refer  to  specific  one-electron 
spin  orbitals  (with  |  i,j,k ...  j  and  { a,b,c, . . .  ( 
labeling  occupied  and  virtual  orbitals,  respectively). 
Then  the  matrix  element  is  represented  convention¬ 
ally  by 


.  t 

(*r  1 1.2 . In  —2),a,j  )  ^-a/tl,2 . in -2),i,b  , 


where  orbital  a(b)  has  replaced  orbital  Hj).  Since 
iJ+\  1,2. . .  (n  —2),a,b  ), 

•  (*^)— v(;|«')  ’ 146 

+i,>  (''«*)  ■ 


The  sign  of  the  first  term  arises  from  the  permuta¬ 
tion  of  orbitals  required  to  match  up  the  remaining 
orbitals  when  /  =?=/.  When  i  no  rearrangement 
is  necessary  (second  term).  Again,  the  diagonal 
matrix  elements  are  zero  if  $  is  real.  In  second 
quantization,  we  respresent  3/3 R  as  a  one-electron 
operator: 


■  <47> 

where 


Then,  following  the  usual  procedures  (and  indicat¬ 
ing  only  the  orbitals  which  differ  from  those  in 

<t>„> 

(fa  ~dRb'J'} 

=  &ij<}ob  —Sab4ji  + 

s 

-8»(“!m‘)-5-‘(j jw')- 

(49) 

The  nonzero  configuration  derivatives  within  the 
chosen  model  space  are  shown  digrammatically  in 
Fig.  1.  In  these  diagrams  an  open  dot  represents 
the  differential  operator,  which  (because  of  the  or¬ 
bital  product  form  of  the  model  functions)  acts  as 
one-electron  operator.  Thus,  each  open  dot  must 
have  one  vector  (a  particle  or  hole  line)  entering 
and  one  leaving.  Particle  (hole)  lines  are  oriented 
from  right  to  left  (left  to  right).  The  terms  re¬ 
quired  are  specified  diagrammatically  by  forming 
all  possible  connections  of  the  bra  and  ket  model- 
state  vertices  with  the  differential  operator.  The 
bras  and  kets  are  indicated  by  diamonds  with  equal 
numbers  of  specifically  labeled  particle  and  hole 
lines  (the  so-called  “open”  lines)  which  describe  the 
model  states  with  reference  to  the  Hartree-Fock 
determinant.  The  algebraic  expressions  are  re¬ 


gained  diagrammatically  by  interpreting  the  deriva¬ 
tive  interaction  symbol  as  (outgoing  vector  1 3/3 R 
incoming  vector)  with  a  sign  factor  ( —  1  )/+*, 
where  /  indicates  the  number  of  loops  and  h  the 
number  of  hole  lines.  Thus  the  two  terms  of  Eq. 
(46)  are  represented  by  the  SE  j  SE  diagram  of  Fig. 

1  and  its  particle-hole  reversed  (phr)  diagram. 

Since  the  integrals  in  Eq.  (46)  involve  two  orbitals 
of  the  same  type  (occupied  or  virtual),  the  entering 
and  leaving  vectors  of  the  differential  operator  go 
in  the  same  direction. 

The  derivative  parts  of  the  j?'"1  terms  in  Eq. 

(40)  are  diagrammed  and  derived  similarly.  The 
usual  two-electron  interactive  perturbation  V 
[represented  by  a  solid  dot  and  referred  to  as  Vs  in 
Ref.  (21)]  in  the  right-hand  pan  of  the  first  term 
must  have  two  vectors  entering  and  two  leaving. 
Figure  2  shows  all  of  the  types  of  diagrams  which 
occur  in  the  evaluation  of  this  term  for  a  model 
space  consisting  of  the  Hartree-Fock  determinant 
plus  selected  singly  and  double  excited  deter¬ 
minants.  The  diagrams  for  the  second  term  are 
obtained  by  time  reversal  (TR)  (left-right  mirror 
imaging)  of  the  ones  shown.  The  off-diagonal 
1  block  diagrams  will  come  from  TR  of  their  conju¬ 
gate  blocks;  e.g.,  TR(SE|DE)  yields  the  DE  SE  di¬ 
agrams.  We  note  that  the  DEISE  terms  in  the  fig¬ 
ure  can  have  triply  excited  intermediate  states 
i  (whereas  the  SE|DE  terms  cannot)  because  of  the 
ordering  of  the  two  operators.  Of  course  the  oppo- 
i  site  holds  when  the  order  is  reversed. 

!  The  diagrams  in  Fig.  2  may  be  derived  algebrai¬ 
cally  in  the  following  manner.  Let  us  examine  the 
|  case  of  a  SE|SE  (single-excitation  diagonal  block) 
matrix  element.  Take  the  definitions  in  Eq.  (44) 
and  let  <PM  be  a  single  excitation 

|  <t>M)=ctk  |  <t>o>  •  '  (50) 


Then 


FIG.  1.  Zeroth-order  diagrams  for  the  first  derivative 
part  of  the  nonadiabatic  coupling  matrix  elements. 


I  (51) 

Next,  we  determine  that 

<<t>„!  V*t)m-(bj\  \ck)  ,  (52) 

since  6*  SJk  is  not  allowed  (i.e.,  <t>M  must  not  be  a 

model  state).  The  combination  of  Eqs.  (51)  and 

(52)  yields  two  types  of  terms  ,  i 
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FIG.  2.  First-order  diagrams  for  the  first-derivative  part  jfj  of  the  nonadiabatic  coupling  matrix  elements. 


4»r  <  <I>M  I  ^4»5>  =  2  ^  I  |K>  . 


which  are  graphically  given  by  the  first  SEjSE  diagram  of  Fig.  2  and  its  phr.  The  remaining  orbital  sum¬ 
mation  corresponds  in  the  diagram  to  the  vector  not  touching  a  model-state  vertex  (indicated  by  a  diamond). 
Similarly,  if  <t>M  is  a  double  excitation 

|  «t>(,>*ctd+/Ic  |  <J>0)  »  (54) 

then 

a**")  =  ^  (y11’2 . . (»  —  3),i,c,d  (55) 

may  be  reduced  to 


(56) 


where  the  notation 


is  used.  Taken  together  with 

<4>r| F<t>s>=5 >j«iAlk)b  |  \cd)+8blcd)((dc)j\  \kl >  , 
Eq.  (56)  gives 


(*r  I  K<t>s)  =  _S>/2  (™  i  \oe)+&ab^,^n  |  in 

+  2  (j  ^-^<'*1  i,m>  • 


Y  jaFy</i>l  ^ae'>~'^‘\n\i'bRb/  <a;l  '  (59) 

The  first  two  terms  of  Eq.  (59)  correspond  to  the  second  diagram  (in  Fig.  2)  and  it  phr.  The  last  two  are 
represented  by  the  remaining  SEJSE  digram  and  its  phr. 

The  second-derivative  matrix  elements  may  be  derived  similarly.  From  the  form  in  Eq.  (43)  we  substitute 
for  Glnl  and  obtain  at  low  order 


WiM-gr*})- 
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The  complexity  of  successive  terms  grows  rapidly;  the  second-order  part  consist  of  17  terms. 

Let  us  convert  the  zeroth-order  contribution  to  orbital  form.  In  second  quantization  we  represent  3:/3 R2 
by  3/3 R  operating  twice: 

Then  we  find  for  single  excitations,  for  example,  that 

+2S“f(’|af')(y  as')+5,*([af'  las') 


These  terms  are  prescribed  diagrammatically  by  the  SE|SE  block  of  Fig.  3.  Diagrammatically,  the  necessary 
terms  consist  simply  of  all  possible  proper  connections  of  the  SE  bra  and  ket  with  two  derivative  operator 
symbols.  Note  that  the  intermediate  states  here  can  be  model-space  states  since  they  are  part  of  the  resolu¬ 
tion  of  the  identity  in  the  form 

l=*2i(><fl  •  <63) 


The  factors  of  2  on  the  last  two  terms  of  Eq.  (62)  arise  diagrammatically  when  it  is  observed  that  the  last 
two  diagrams  are  equivalent,  as  are  the  third  from  the  last  and  its  phr.  The  remaining  diagrams  are  only 
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reconcilable  with  Eq.  (62)  after  combining  the  terms  with  matching  delta  functions.  This  is  accomplished 
by  inserting  the  resolution  of  the  identity  into  the  even  numbered  terms,  and  separating  the  sums  over  occu¬ 
pied  and  virtual  orbitals;  e.g., 


The  first  two  terms  cancel  leaving 


which  then  corresponds  exactly  to  the  first  SE|SE 
diagram.  The  other  two  diagrams  are  obtained 
similarly.  Comparison  with  a  conventional  deriva¬ 
tion  is  straightforward,  and  uses  the  anti- 
Hermitian  property  of  3/3 R;  e.g.. 


The  elementary  second  derivative  terms 


which  are  similar  to  the  £',0)  matrix  elements,  are  j . 
the  only  additional  component  necessary  to  evalu-  | ! 


ate  the  first-order  terms.  The  integral  and  en¬ 
ergy  (first)  derivatives  are  the  same  ones  already 
used  in  E  ,s.  (34)  and  (41),  but  this  time  they  are 
multiplied  by  £'l0'-like  terms.  These  terms  are 
the  ones  required  for  £,u’  in  Eq.  (40).  £''2)  will 
require  second  derivatives  of  the  two-electron  in-  i 
tegrals.  but  these  are  also  available.12 

If  the  nuclear  Schrodinger  equation  is  to  be 
solved  exactly,  either  (1)  the  couplings  must  be  cal¬ 
culated  in  the  coordinate  system  in  which  the 
scattering  problem  is  to  be  solved,  which  is  not  the 
usual  case,  or  (2)  a  complementary  set  of  coupling 
matrix  elements  of  similar  one-electron  operators 
must  be  provided  to  enable  a  transformation  of 
coordinates.3  The  diagrams  presented  in  this  paper 
are  sufficient  in  either  case.  As  derivatives  with 
respect  to  various  coordinates  are  called  for,  one- 
electron  equations  analogous  to  Eq.  (47)  may  be 
defined.  However,  the  diagrams  required  will  be 
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FIG.  3.  Zeroth-order  diagrams  for  the  second  derivative  pan  Q"  of  the  nonadiabatic  coupling  matrix  elements. 
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the  same  for  any  one-electron  operator.  A  comple¬ 
mentary  set  of  couplings  may  be  obtained  in  the 
same  fashion.  Similarly,  cross-term  matrix  ele¬ 
ments  arising  in  the  transformation  can  be  evaluat¬ 
ed  using  the  second-derivative  diagrams  if  the  two 
one-electron  interactions  are  appropriately  labeled. 
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APPENDIX:  AN  EXACT  MATRIX 
HELLMANN-FEYNMAN  RELATION  | 


Write  the  eigenvalue  problem  in  the  chosen 


model  space 

HM=M 

(Al) 

(with  £  diagonal)  as 

b'hb=E  . 

(A2) 

( HB  can  be  H£_  of  from  the  text.)  Then  dif¬ 


ferentiate  each  side  of  Eq.  (A2)  with  respect  to  R 
'representing  this  by  a  prime  symbol): 


(Mi)'M+Bifi'i+B*M'=£'  • 

(A3) 

Using  Eq.  (Al)  and  reorganizing,  we  obtain 

B*U'B  =£■-(£')■  BE -ZB* i'  . 

(A4) 

From  differentiating  the  normalization  equation 

5^=1, 

(AS) 

we  can  see  that 

(B*VB=-ifZ’  ■ 

(A6) 

Realizing  that  £  and  thus  £’  are  diagonal,  we  ob¬ 
tain  from  Eq.  (A6)  and  the  off-diagonal  part  of 

Eq.  (A4) 

[B'H’BaB\ 

14  4  U- 

(A7) 

When  B  is  real. 

Bay—Bya  . 

(A8) 

Then  using  Eq.  (A8)  in  Eq.  (A6)  gives 

^ayBya— ~^Bya(BarV  , 

r  r 

(A9) 

so  that 

1  M'B'  )ao=0  . 

(A10) 
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